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Abstract 

The static potential in the confinement "phase" of the SU(2) Higgs model is 
studied. In particular, the observation of the screening (called string breaking) 
of the static quarks by the dynamical light quarks leading to the formation of 
two static-light mesons was not observed before my work in non-Abelian gauge 
theories. The tool that I employ is lattice gauge simulation. The observable 
from which the spectrum of the Hamiltonian in presence of two static quarks 
can be extracted, is a matrix correlation whose elements are constructed not 
only from string-type states represented by Wilson loops (like in pure gauge 
theories). Additional matrix elements representing transitions from string-type 
to meson-type states and the propagation of meson-type states are taken into 
account. From this basis of states it is possible to extract the ground state and 
first excited state static potentials employing a variational method. The crossing 
of these two energy levels in the string breaking region is clearly visible and the 
inadequacy of the Wilson loops alone can be demonstrated. I also address the 
question of the lattice artifacts. For this purpose lines of constant physics in the 
confinement "phase" of the model have to be constructed. This problem has only 
partially been solved. Nevertheless it is possible to show that the static potentials 
have remarkable scaling properties under a variation of the lattice spacing by a 
factor two and are almost independent of the quartic Higgs coupling. 
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Chapter 1 
Introduction 



The success of the quark- constituent picture both for resonances and for deep- 
inelastic electron and neutrino processes makes it difficult to believe quarks do 
not exist. The problem is that quarks have not been seen. 



K.G. Wilson, 1974 



The theory of strong interactions plays a pivotal role in particle physics. It is 
part of the Standard Model of elementary particles which successfully describes 
the constituents of the matter in terms of quantum gauge field theories. These 
theories are based on the gauge principle 0] : the fields in the theory have internal 
degrees of freedom associated with a gauge group and it is required that local 
transformations of these degrees of freedom leave the physics unchanged. The 
gauge group of the Standard Model is SU(3)^ x SU(2)l x U(1)y and the degrees 
of freedom associated with them are color for SU(3), weak (left-handed) isospin 
for SU(2) and hypercharge for U(l). The gauge group together with the gauge 
principle dictate the structure and properties of the interactions. The particle 
content, described by means of relativistic local quantum fields, has to be deduced 
from what nature tells us. 

The particles which take part in strong interactions are called hadrons: Gell- 
Mann and Zweig 0, |^ proposed a model that explained the low energy prop- 
erties of the hadrons (like mass and spin) in terms of elementary constituents 
called quarks. Bjorken [Q studied, within the framework of current algebra, the 
electron-nucleon scattering and discovered the scaling property of the structure 
functions for large electron momentum transfer (deep inelastic scattering). The 
Bjorken scaling was experimentally confirmed and could be understood with the 
assumption that the electrons scatter off almost-free pointlike constituents 
inside the nucleon, which were called partons |^, 0. Later the partons were iden- 
tified with the quarks on the basis of their quantum numbers. The question at 
that moment was to find a theory in which particles are free at high energies. The 



decisive step was then made with the proof that non-Abehan gauge field theories 
exhibit asymptotic freedom p|, ^, 0]. The strength of the interaction given by 
the gauge couphng becomes weak at shorter distances (or equivalently at high 
energies) and this is consistent with the Bjorken scahng. In order to resolve 
several difficulties of the quark model, like the construction of an antisymmetric 
wave function for the A"*""*" baryon and the discrepancy between the prediction 
and experimental data on the total cross section e"*"e~ —>■ hadrons, it was al- 
ready suggested that quarks must have a new quantum number called color and 
exhibit color symmetry. Fritzsch, Gell-Mann and Leutwyler |ll|] proposed that 



the theory describing the quark dynamics is a non-Abelian gauge theory with 
gauge group SU(3) associated with the color symmetry. This theory was named 
quantum, chrom,odynam,ics (QCD): its ingredients are quarks and gluons, usually 
called partons. The gluons are the vector bosons that mediate the interactions: 
in contrast with an Abelian gauge field theory where the vector boson (the pho- 
ton) is gauge neutral, the gluons carry color quantum numbers and therefore have 
self-interactions. It is this property which is responsible for asymptotic freedom. 
Due to asymptotic freedom, the short distance behavior of the partons can be 
described with a perturbative expansion in the small value of the gauge coupling. 
Within the perturbative approach, QCD found important confirmations as the 
theory of strong interactions, such as the prediction of a logarithmic deviation 
from Bjorken scaling in structure functions, confirmed experimentally in deep 
inelastic lepton-nucleon scatterings. 

What is observed in nature are not the partons, but the hadrons, which 
are color-neutral objects. The fact that colored partons cannot be seen iso- 
lated led to the conjecture of color confinement the partons are always bound 
into the hadrons. In order to prove this assumption from QCD one should be 
able to describe its properties at long distances corresponding to the size of the 
hadrons. Perturbation theory is not applicable because the gauge coupling is 
large at this scale. Wilson |T^ proposed in 1974 a new approach to gauge field 
theories amounting to the discretisation of the four- dimensional space-time on a 
Euclidean lattice. The quantisation of this theory is naturally performed in the 
path integral formalism. The matter fields are treated as classical variables liv- 
ing on the points of the lattice and the gauge field is represented by connections 
(links) between the matter fields on nearest-neighbor points. The quantum effects 
in the observables of the theory are introduced by evaluating their expectation 
values expressed as Feynman path integrals |T^. In the Euclidean lattice formu- 



lation, a quantum field theory looks like a classical statistical system. Particle 
and solid state physics mutually profited by this relationship [Q . The concept of 
renormalisation of a gauge field theory receives new insights. The regular isat ion 
of the theory on the lattice is associated with an ultraviolet cut-off, the inverse 
lattice spacing a~^. The field theory is changed in the short distance region while 
its long distance properties are preserved. The question one is interested in, is 
whether there is the possibility of constructing a continuum quantum field theory 



c 





Figure 1.1: The Wilson loop, here represented on a lattice, describes the confine- 
ment of static quarks. 



from the lattice field theory: that is, is the limit a — i> of the lattice field theories 
well defined? To answer this question, we should be able to reproduce the same 
physical situation on lattices with different cut-offs a and consider the behavior 
of dimensionless physical quantities when a — *> 0. The equations describing the 
change in the parameters of the theory under variation of the lattice spacing a 
are the renormalisation group equations (RGBs) of the lattice field theory. One 
consequence of the RGEs is that the continuum limit of a lattice regularised 
asymptotically-free theory is reached when the lattice bare gauge coupling g is 
sent to zero. To show the correspondence between a lattice field theory and a 
statistical system one considers the field propagator on the lattice. For example, 
in a statistical system of Ising spins, the corresponding quantity is the spin-spin 
correlation, whose exponential decay is governed by the correlation length. On 
the lattice, the correlation length equals the inverse mass gap. By keeping the 
mass gap fixed at its physical value, the correlation length expressed in units of a 
diverges in the continuum limit. Thus, the continuum limit of a lattice field the- 
ory, if it exists, corresponds to a second order phase transition in the parameter 
space of the statistical system. 



Wilson |T2[ originally proposed lattice gauge theories in order to explain color 
confinement. To this end, he derived an expansion valid for strong gauge coupling 
in which confinement arises naturally. However, in non-Abelian gauge theories 
the continuum limit is reached when g ^ due to asymptotic freedom. Another 
method must be developed to study the confinement in the weak gauge coupling 
regime. We consider the system composed of a pair of infinitely heavy or static 
quark and anti-quark. The static quark (anti-quark) is treated as an external 
source in the (complex conjugate of the) fundamental representation of the gauge 
group. In pure SU(A^) Yang-Mills gauge theories, the potential between the static 
quarks, called the static potential, can be extracted in the path integral formalism 
from the expectation value of Wilson loops represented in Fig. |1 . 1| . On the lattice, 
they are defined as the trace of the product of the gauge hnks U over a closed 



path composed of two straight tiine-hke hnes and arbitrary space-hke paths C 
and C connecting the time-hke hnes: 

(l^cc) = (trH^) ~ e-'-^QQ^^-y) ^ (1.1) 

where t is the time extension of the Wilson loop and VQQ(r) is the static potential 
for the separation r of the static quarks. The expectation value in eq. ( p..l|) can be 
computed by Monte Carlo simulation of the Yang-Mills theory on the lattice. The 



seminal work was done by Creutz |jT5[ for the gauge group SU(2) and since then, 
there have been a number of detailed studies which show a linear confinement 
potential dX large distances between the static quarks close to the continuum limit, 
both for gauge group SU(2) [|l|, and SU(3) P, |19|. 



When the Yang-Mills gauge theories are coupled to matter fields in the fun- 
damental (quark) representation of the gauge group, the potential between a pair 
of static quarks is expected to flatten at large distances: the ground state of the 
system is better interpreted in terms of two weakly interacting static-light mesons 
which are bound states of a static and a dynamical quark. The dynamical quarks 
are pair-created in the strong gauge fleld binding the static quarks. This phe- 
nomenon is called string breaking or screening of the static charges. The name 
"string" refers to the gauge fleld conflguration which conflnes the static quarks 
and leads to the linear conflnement in pure gauge theories. 

In recent attempts in QCD with two flavors of dynamical quarks, this string 
breaking effect was not visible |^, ^, |2^, ^, ^ |2^. The string breaking distance 



Tb around which the static potential should start flattening off, could nevertheless 
be estimated in the so called quenchedQ approximation of QCD to be 



r-b ~ 2.7ro, (1.2) 



where the scale tq ~ 0.5 fni was introduced in |17|. The static potential at short 



distances and the mass of the static-light meson can be computed in quenched 



QCD. The approximate value in eq. ( |1.2| ) was obtained from the crossing point 



of the linearly rising potential with twice the value of the meson mass (which is 
expected to be the asymptotic value of the potential after string breaking). 

The investigation of the static potential in models other than QCD is therefore 
relevant in order to understand its origin and identify possible failures of the 
methods used to extract it. First studies of string breaking were performed 
with a hopping-parameter expansion in SU(2) gauge theory with Wilson fermions 



2^. In the Schwinger model, which is quantum electrodynamics (QED) in two 
dimensions, the exact solution for the static potential can be given in the limit of 
zero fermion mass [E3]: V{r) = {e^/^T / 2) {1 —exp{—er / ^/n)} , where e is the charge 



of the static sources. String breaking was established by numerical simulation in 



^ In this approximation the effects of internal quark loops are neglected. In practical Monte 
Carlo simulations the computational effort is considerably reduced. 



the Schwinger model ^^ |30[ . Numerical evidence of the screening of the static 
potential was also found in the U(l) Higgs model (scalar QED) in two dimensions 
^]. The flattening of the static potential at large distances is also expected in 



the confinement "phase" of the SU(2) Higgs model. Indeed, early simulations 



yielded some qualitative evidence for string breaking ^, ^ . 

String breaking can also be studied in Yang- Mills theories using static sources 
in the adjoint representation of the gauge group. The gauge field itself is re- 
sponsible for the screening of the sources and the formation of hadrons called 
"gluelumps". An important numerical investigation concerning this screening 



has been carried out by C Michael in |34], where it has been noted that string 



breaking can be a mixing phenomenon. The static potential is extracted from a 
matrix correlation in which two types of states enter, the adjoint "string" and the 
"two-gluelump" . However, due to large errors, no clear evidence for string break- 
ing could be given. The first numerical evidence using the mixing method for 
string breaking in non-Abelian gauge theories with dynamical matter fields, was 
given in the four-dimensional E^ and three-dimensional PBI SU(2) Higgs model 



by the computation of the potential between static quarks. Most recently, the ex- 
traction of the static adjoint potential in the three-dimensional SU(2) Yang- Mills 
theory |3^, |38| shows also evidence for string breaking. 

Finally, we want to mention that string breaking has been seen in finite tem- 
perature QCD [^, where the static potential can be extracted from Polyakov 
loop correlators. 



The status quo before our work, was that no clear evidence for string breaking 
in non-Abelian gauge theories was established. In our research, we investigate 
the potential between static quarks in the four dimensional SU(2) Higgs model 
on the lattice. In Chapt. |^, we describe the model. The parameter space of the 
theory is divided in two "phases" , the confinement and the Higgs "phase" . In the 
confinement "phase", the properties are similar to QCD: screening of external 
charges by the dynamical Higgs field is expected. We describe the error analysis 
of the statistical measurements. 

In Chapt. 1^, we concentrate on the determination of the mass spectrum of the 
static-light mesons, which are expected to be the asymptotic states after string 
breaking. We describe the variational method that we use for extracting the 
energy spectrum from a matrix correlation function constructed with a basis of 
states that can mix. We will use the same method for the determination of the 
static potential. The basis of states is enlarged by the use of smeared fields: we 
present a study of different smearing procedures for the Higgs field. 

In Chapt. ^ we introduce the matrix correlation function from which we ex- 
tract the static potential. We use two "types" of states: "string states" and 
"two-meson states". The variational method determines the best linear combi- 



nation approximating the true eigenstates of the Hamiltonian. We present the 
results for P = A/ g^ = 2.4, where g is the gauge couphng. We are able to deter- 
mine the ground state and first excited state static potential with good accuracy. 
This allows us to study the overlaps of the approximate eigenstates, determined 
from our basis of states, with the true eigenstates. For comparison with the re- 
cent studies of string breaking in QCD, we analyse what happens if we use only 
the Wilson loops for determining the ground state static potential. 

In Chapt. ^, we address the question of the "continuum limit": in order to 
investigate lattice artifacts in our results, we would like to reproduce the physical 
situation of Chapt. ^ on a coarser lattice at /? = 2.2. In the parameter space 
of the model, this would define a line of constant physics (LCP), along which 
two dimensionless physical quantities are kept constant under variation of the 
lattice spacing. The static potential provides us with a first quantity sensitive 
to the mass of the dynamical Higgs field. We study the definition of a second 
quantity sensitive to the quartic Higgs coupling. Although we are not able to 
match precisely the parameters along the LCP, we find a parameter region in 
which the discussion of the scaling properties of different quantities, in particular 
the static potentials, is possible. 

In Chapt. 1^, we summarise the results of our work and give some prospectives 
for future investigations. 

A number of more technical information is relegated to appendices. In Ap- 
pendix 0, we explain the notation conventions that we use throughout the work. 
In Appendix ^, we construct the transfer matrix operator for the SU(2) Higgs 
model and prove its positivity, which is the condition for a real energy spectrum 
of the theory. The connection between path integral expectation values and vac- 
uum expectation values of corresponding time ordered operators is also shown. 
In Appendix P, we describe our algorithms for the Monte Carlo simulation of 
the SU(2) Higgs model. In Appendix 0, we explain the implementation of the 
one-link integral method which allows the reduction of the statistical variance of 
the correlation functions. Finally, Appendix ^ is devoted to the description of 
the parallelised computer program that we use for the Monte Carlo simulations. 



Chapter 2 

The SU(2) Higgs model 



2.1 Definition of the model 

The SU(2) Higgs model on a four-dimensional Euclidean lattice^ is defined by 
means of a gauge field of dimensionless SU(2) matrices f/(x, /x) and a complex 
Higgs field <l>(x) in the fundamental representation of the gauge group SU(2) and 
with canonical mass dimension one. The full action is 



T. ^11. 



+ ml^''{x)^{x) + Xq $"I'(2;)$(x) 



(2.1^ 



where Sw is the Wilson action eq. ( |A.8D for the SU(2) gauge field. Introducing 
the dimensionless lattice fields $l = {(^/V^)^ (we drop the subscript L in the 
following) together with the new couplings A = k^Xq and k = (1 — 2X) / {8 + a^ml) 
the action can be written as 



(2.2) 



S = 5w + E j '^\x)<!>{x) + X [$t(x)$(x) - 1 



The physics of the model is controlled by the three dimensionless bare parameters 
(3 = 4:/g'^, K,, X. We will use the parametrisation eq. ( ^.2| ) throughout the work. 
We can rewrite eq. ( p.2| ) using the 2x2 matrix notation for the Higgs field 



(p{x) = p{x)a{x), p{x) > 0, a{x) G SU(2) defined in Sect. |A.3 



5 = ^w + E {^ tr {v\xMx)) + A [i tr {^H^Mx)) - 1 
— ^E t^ i^H^W (x , fi)(p{x + ajl)) > . 



(2.3) 



^ For a detailed description of the notation we refer to appendix ^. 
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The gauge and Higgs field are both represented as 2 x 2 matrices which are 
equal to a real constant times an SU(2) matrix. This is particularly useful for 
programming purposes, see Sect. [E.2| . 

The Euclidean expectation value of an observable^ 0\U, $] is defined by the 
Feynman path integral 

(O) = i|D[t/]D[0]O[f/,$]e-^[^'*], (2.4) 

where Z is the partition function 

Z = I V>[U]T)[<P\e-^^^^''\ (2.5) 

D[f/] denotes the product measure df/(x, jj) (df/ is the Haar measure on SU(2)) 
and D[0] denotes 11 d0i(x) ■ ■ ■ d04(x). The fields 0i(x) {i = 1, 2, 3, 4) are the four 

X 

real components of ^(x) defined in eq. ([A.14|) . 



The expectation values in eq. (|2.4| ) can be shown to correspond to Euclidean 
vacuum expectation values of corresponding time-ordered operators for large 
enough time extent of the lattice. This correspondence can be established through 
the definition of a time evolution operator called transfer matrix. How this is done 
for the SU(2) Higgs model is the subject of Appendix |B[ 

The lattice formulation of a quantum field theory provides a mathematically 
well-defined, non-perturbative and completely finite regularisation of the theory. 
An analytical solution of eq. ( p. 41) is in the most cases not possible, but it can 



be computed with Monte Carlo algorithms. The overall aim |40] is to generate 
a representative ensemble of field configurations {[t/„,<l'„], n = 1,2,...,A^} for 
the path integral eq. ( |2.4| ) by employing a stochastic process. Representative 
means that the probability distribution of the configurations in the ensemble is 
the Boltzmann distribution exp(— ^[f/, $]). We evaluate then eq. ( |2.4| ) as the 
ensemble average 

_ 1 ^ 

O = T7E0[f/n,$n]. (2.6) 

The value O has a statistical error A(0). Moreover, in most cases one is interested 
in secondary quantities, which are functions of the primary averages O. 

The methods that we use for the Monte Carlo simulation of the model are 
described in detail in Appendix |C|. The determination of the statistical errors 
will be the subject of Sect. 



^ We put square brackets to denote the dependence on the whole field configuration and not 
only on a particular field variable. 



2.2 Symmetries 

We dedicate a section to the discussion of symmetries of the action and integration 
measure which reflect themselves in useful properties of the expectation values 
defined by eq. (^). 

Let us start from the four-component 0^ theory without gauge interactions. 
The four components of the scalar field 4>i{x) {i = 1, 2, 3, 4) can be put in a 2 x 2 



matrix as described in Sect. |A.3| . In this notation the action of the pure scalar 
theory is 

M M -|2 

-K^ tr(v9"^(x)v?(x + a/i))|. (2.7) 



The model is symmetric under the global 0(4) transformation (f)i{x) = Rij(j)j{x), 
where i? G 0(4). In the matrix notation for the scalar field this transformation 
is equivalent to a global 811(2)^ ® SU(2)j^ symmetry defined as 

ifix) -^ Ay{x)A^ , Al,rGSU(2). (2.8) 

The SU(2) Higgs model is obtained by gauging the subgroup 811(2)^, i.e. by 
"promoting" the global symmetry associated with the group 811(2)^ to a local 
symmetry. According to the gauge principle ||l|, this automatically requires the 
appearance of the gauge field U{x, fi) and of interactions. We end up precisely 
with the action in eq. ( |2.3| ). Under gauge transformation defined by the field of 
SU(2) matrices {A(x) = Al(x)} the fields transform as 

U^{x,n) = A\x)U{x,fi)A{x + afi) (2.9) 

if^ix) = A\x)ip{x). (2.10) 

The SU(2) Higgs model has a residual global SU(2) symmetry defined by the 
diagonal subgroup Aj^ = Ar = A: 

ip{x) — > A^ip{x)A, (2.11) 

f/(x,/i) — > A^U{x,iJ.)A. (2.12) 

This symmetry is called the (weak) isospin. 

Let us now discuss the symmetry under gauge transformation. The action 



of the SU(2) Higgs model eq. ( p.3|) is invariant under the gauge transformations 
eq. ( |2.9|) and eq. ( p.lO|) per construction. The integration measure in eq. ( |2.4|) as 
well: for the gauge field the property dU = dU^ is a consequence of the invariance 
of the Haar measure: dU = d{UV) = d{VU) for V E SU{2). For the Higgs field 
we note that 

/ d0i(x) ■■ ■d04(a;) oc / p^dp da, (2.13) 

JR4 Jo JSV{2) 



if we write the Higgs field in the notation of eq. ( [A.15| ) (f{x) = p{x)a{x). From 
(p^{x) = p{x){A^ {x)a{x)) the invariance of the measure follows immediately. 
Defining 0^[U, $] = 0[U^, $^] as the gauge transform of the observable O, we 
therefore have the property 

(O^) = (O). (2.14) 

In the continuum SU(2) Higgs model one speaks of Higgs and W-boson fields 
which are not gauge invariant: a particular gauge is fixed and perturbation theory 
can then be applied. In the non-perturbative lattice scheme one does not need to 
fix the gauge and gauge-invariant definitions of interpolating fields for the Higgs 
and W-bosons have then to be used. A gauge invariant composite Higgs field is 
defined as 

H{x) = ^\x)^{x) = -tT{ip\x)ip{x)). (2.15) 

It is an isospin scalar field. A gauge invariant W-boson field is defined as 

Wr^ix) = -ttT{TrV^{x)) (r = 1,2,3), (2.16) 

Vfj_{x) = a^(x)f/(x,/i)a(x + a/i) , (2-17) 

where V^(x) is called the gauge invariant link variable. The field in eq. (|2.16| ) is 
an isospin 1 (with isospin index r) spin 1 vector field. The isospin property can 



be seen transforming the field under the isospin transformations eq. (|2.11|) and 



eq. ( |2.12| ) and using the relation 

ATrA^ = R{A)srTs, A G SU(2), i?(A) G S0(3) , (2.18) 

which defines R{A) as the adjoint (isospin=l) representation of the isospin group 
SU(2). 

Another symmetry of the action and the integration measure is under com- 
plex conjugation of the field variables, Uab{x) — >■ t/*f,(x), ^a{x) -^ ^l{x). For 
the measure this symmetry is a consequence of the equivalence of the repre- 
sentation 2 and 2* of SU(2) eq. ( A.13| ). For an observable of the form L 



^^{x)U{x, y)^{y), where U{x, y) is a link path connecting y with x, we can write 
L* = {^*{x)yU*{x,y)^*{y) from which it follows 

in = (L). (2.19) 

The reality of the action and the measure gives the result (L) G R. 

2.3 Phase diagram 

The bare parameters k, A are restricted to the ranges k > 0, A > 0. A negative 
value of A would lead to an action which is not bounded from below and therefore 
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K 



pure gauge 

(K=0) 



Higgs "phase" 
1st order p. t. 



confinement "phase" 
( ~ QCD ) 



2nd order p.t. 



(j)"^ theory 
(g = 0) 



Figure 2.1: The phase diagram of the SU(2) Higgs model in the (/3, k) plane is 
shown for a typical value of A. 



the path integral eq. 
transformation 

$(a;) 



41) would not be well defined. With the help of the 



-$(x) if 
^(x) if 



.^Em^''/'' = _i 



(2.20) 



it is easy to see that the partition function eq. ( |2.5| ) has the following symmetry 

Z{f3,K,X) = Z{f3,-K,X). (2.21) 

Replacing k with —k, corresponds merely to a change in the description of the 
model and leaves the physics of the model invariant: with the help of the trans- 
formation eq. ( |2.20| ) we can define a mapping of the observables of the model such 
that the expectation values eq. ( p. 41) computed with parameter k are reproduced 
by expectation values computed with parameter —n. Without loss of generality 
we can therefore restrict ourselves to the parameter range k, > 0. 

IA\ shows the phase structure of the SU(2) Higgs model |^ ^ 



Fig. 



In the 



{/3, k) plane there is a line which is believed to be a first order phase transition line 
separating the confinement "phase" at small values of k, from the Higgs "phase" 
at larger k. The situation in Fig. 2J. is for a typical value of A. The /? = oo 
{g = 0) boundary of the phase diagram is the pure scalar 0^ model, in which there 
is a second order phase transition line Kc{X) separating the phase where the 0(4) 
symmetry is spontaneously broken [k > Kc) from the symmetric phase [k < Hc)- 
At K = the Higgs field becomes infinitely heavy and decouples from the gauge 
field: we are left with a pure SU(2) gauge theory. At small values of /3 there is 
an analytic connection between the two "phases" : this is why we put "phase" in 
quotation marks. As an analogy, we can think of the liquid and vapor regions in 
the phase diagram of a fiuid. Nevertheless, the physical properties can be quite 
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different in the two "phases" , as can be seen for instance by inspecting the static 
potential extracted from the Wilson loops at small distances. In the confinement 
"phase" , which is the continuation at finite /3 of the symmetric phase of the (f)'^ 
theory, the gauge field has confinement properties like in QCD. The static charges 
are bound by the color field string and the potential rises, in good approximation, 
linearly as in pure gauge theory ||32|, ^. In the Higgs "phase", which continues 
the spontaneously broken phase of the 0^ theory, the Higgs mechanism is at 
work: the gauge vector W-bosons become massive. Far enough from the phase 
transition the interaction between static charges is mediated by the exchange of 



W-bosons and the potential has a Yukawa form [42, 41 



Because we are interested in the Higgs model as a test model for QCD, we 
work in the confinement "phase" and at large values of the gauge coupling g 
{P = A/g^ = 2.2, 2.4). The situation of the Standard Model Higgs sector is 
"opposite", in the sense that it is in the Higgs "phase" and natural choices for 
the gauge coupling are /3 ~ 8 [^. The simulation in the confinement phase 
at these small values of the gauge coupling would be impossible because of the 
extremely large correlation lengths. Close to the continuum, in the lowest order 
approximation of the massless perturbative renormalisation group equation, the 



lattice spacing a depends exponentially on P [41 



In (cAl) 



12 
43 



vr^/?: 



(2.22) 



where Al is the renormalisation group invariant A-parameter in the lattice scheme. 
The scale of a lattice gauge theory simulation can be set by computing the phys- 
ical length To [[1^ from the force between static charges. In QCD this length 
corresponds to 0.5 fm. In Sect. [4.2| , we will present the results of the simulation 
of the SU(2) Higgs model for the parameter set (3 = 2.4, k = 0.2759, A = 0.7: 



the scale tq is approximately 5 lattice spacings. If we evolve with eq. (|2.22|) the 
lattice spacing from (3 = 2.4 to /5 = 8 we find 



To 
a 



^ ~ 10^ iP 



(2.23) 



2.4 Monte Carlo simulation 



The general principles of a simulation of a quantum field theory on a space-time 
lattice are explained in reference [Q . Here and in Appendix |C|, we give a detailed 
description of the Monte Carlo updating algorithms that we use for the simulation 
of the SU(2) Higgs model. 

We use a hybrid over-relaxation algorithm (HOR) 



3|, y, ^ |6|, ^ ^ ^, 
SC , ^ ^ which is a mixture of heatbath and over- relaxation algorithms. These 
algorithms are local in the sense that in each step only one field variable ip (a gauge 
link U{x, /i) or a Higgs field variable $(a;)) is updated. A sequence of local steps. 
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updating all field variables, is called a sweep. The updating of the field variables 
is a stochastic process: the change ip' of the field variable ip happens with a given 
transition probability p{%p — »• ip'). In order that the field configurations generated 
reach the equilibrium distribution exp(— S')/Z, where 5* is the action and Z the 
partition function eq. ( |2.5| ), it is sufficient that the updating algorithms satisfy 
the conditions of local detailed balance and ergodicity. Local detailed balance 
means that 

p(^^^')e-^W = p(^' ^ ^) e-^(^') , (2.24) 

where S{ip) is the part of the action depending on the field variable ip. Ergodicity 
means that each field configuration can be reached by a finite number of updating 
sweeps. 

In the heatbath algorithm, the new value ip' for the field variable ip is chosen 
independently of the original value ip according to the transition probability 

In analogy with thermodynamics, we can imagine that the field variable is brought 



in contact with an infinite "heat bath" in the equilibrium distribution eq. (|2.25|) . 
In the over-relaxation algorithm, the new value ip' for the field variable ip 
leaves the action S{ip) invariant (this is called a microcanonical change) 

Sm = SW. (2.26) 

The change ■?/'—>■?/'' is proposed with some arbitrary probability pci'ip —>■ ^') and 
is accepted with probability 

(the factors exp(— 5) cancel if the change is exactly microcanonical). The tran- 
sition probability is p = PaPc- The aim of the over- relaxation is to speed up the 
updating process by choosing the new field variable ip' as far as possible from 
ip (a kind of reflection, see below). A consequence of constant action is that 
the algorithm is non-ergodic. In the HOR algorithms this difficulty is cured by 
combining over-relaxation with heatbath updating sweeps. 

It is often not possible to implement the heatbath and over-relaxation al- 
gorithms exactly. What is then done is to propose a new value ^p' of the field 
variable ip with an approximation of the algorithm and accept the change with a 
probability that corrects for the approximation done. 

The updating of the SU(2) Higgs model we have chosen was inspired by ref- 
erence [^. It consists of cycles, that we call iterations, composed each of one 



heatbath sweep for the gauge field, followed by one heatbath sweep for the Higgs 

13 



field and A'^qr times an over-relaxation block composed by one over-relaxation 
sweep for the gauge field and three over-relaxation sweeps for the Higgs field. The 
integrated autocorrelation times (see Sect. |2.5|) depend on the order in which the 
field variables are updated during the sweep |Q. For the link variables we use 
the SF-updating of |Q, in which the outermost loop runs over the direction /x 
and the internal loops run over the lattice points in lexicographic order (we refer 
to Appendix ^ for the details). The updating sweeps of the Higgs field variables 
process the lattice points in the same lexicographic order. In the HOR algorithms 
the parameter A'^qr should be chosen so as to minimise the autocorrelation times 
of the quantities of interest. Our choice A'or = 1 was motivated by a rough 
study of the integrated autocorrelation times of observables like plaquette, Higgs 
length squared and gauge invariant links. They were found to be minimal for 
A'oR = 1; 2. From the study of these "cheap" (referred to the computer time 
needed for the measurements) observables we could draw useful conclusions for 
the measurements of the observables in which we are interested (see Sect. p.5.2|) . 



In Appendix |y, we give a detailed description of the different parts of the HOR 
algorithm that we use for the simulation of the SU(2) Higgs model. Attention is 
also paid to the generation of random numbers needed for the implementation of 
the algorithms. 

2.5 Statistical error analysis 

An essential part of the Monte Carlo simulations are the estimates of the errors 
of the observables computed as in eq. (f2.6D, called primary quantities, and of 
secondary quantities, which are arbitrary functions of primary quantities. Besides 
the naive statistical error, associated with the finite number of measurements A^ 
and proportional to l/y/N, there are other fundamental sources for errors related 



to the updating algorithms used |55 



• Initialisation bias. The algorithm needs a number of thermalisation steps 
before it "forgets" the arbitrary initial configuration and reaches the ther- 
mal equilibrium where the field configurations are distributed according to 
the Boltzmann factor exp(— 5). 

• Autoccorelation in equilibrium. When thermal equilibrium is reached, the 
field configurations generated by the updating algorithm are correlated. 
This causes the statistical error of O in eq. (|2.6|) to be a factor 2rint(0) 
larger than in an ensemble of independent configurations. The quantity 
'7"int(0) is called the integrated autocorrelation time for the observable O. 

The dependency on the initial (arbitrary) configuration can be avoided by 
waiting a "large enough" number of updating steps before starting the measure- 
ments. By measurements we mean the evaluation of the observables on the field 
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configurations generated by the algorithm. What is "large enough" can be esti- 
mated from the integrated autocorrelation times of the observables. They can be 
very different for different observables. In practice, the observed autocorrelation 
times have almost the same order of magnitude and a number of thermalisation 
steps equal 20 to 100 times the maximum observed autocorrelation time Tint,max 
is a sensible choice. 

The central role in the determination of the statistical errors is played by the 
integrated autocorrelation times. How to estimate them is the subject of this 
section. 

2.5.1 Primary quantities 

We consider a sequence of measurements Ai = A[Ui,^i], i = 1,...,N oi the 
observable A = A[U,^] performed on a large ensemble of field configurations 
{[t/j,$i], i = I, ■■■,N} already in the equilibrium distribution. We denote by 

_ 1 ^ 

^ = l^E^^ (2-28) 

the ensemble average of A. The exact path integral expectation value of A is 
denoted by {A) . If the measurements are statistically independent the value of 
A is normally distributed around the expectation value (A) with variance 



^2 



var(A) = A^-A = {A - A)^ . (2.29) 

The (naive) statistical error is then given by 

(Anaive(A))^ = ^. (2.30) 

In general, there are correlations in the sequence of generated field configurations 
(and hence in the measurements), called autocorrelations and eq. ( p.30| ) under- 



estimates the statistical error. In order to obtain reliable statistical errors we 
follow 0, |§. 

The (unnormalised) autoccorelation function is defined as 

TAit-j) = {iA-{A))iA,-{A)))Mc = TaU-i), (2.31) 

where (■ ■ ■)mc denotes the average over infinitely many independent ensembles of 
configurations in thermal equilibrium. The autocorrelation function F^ depends 
only on the distance between the measurements \t\ = i — j. Typically, it decays 
exponentially 

TAit) ~ exp(-|t|/r) for large t. (2.32) 
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The integrated autocorrelation time is defined as 

-int(^) = ^ E M, (2-33) 

where r^(0) = {{A — (A))"^) = var(y4) is the variance]^ of the observable A. 
Here, "time" refers to the "Monte Carlo time" of the simulation and labels the 
measurements. 

The goal is to estimate the effects of the autocorrelations based on a finite 
(but large) sequence of measurements Ai, i = 1,...,N. The ensemble average A 
in eq. ( |2.28| ) has statistical error, corrected for autocorrelations, given by 

_ 1 ^ 

1 ^"^ / ltl\ 

^^ t=-{N~l) \ ^ V 

^ ^(2rint(A))r^(0) for N :^ t . (2.34) 



Comparing with eq. ( p.30| ), we see that the statistical error is a factor -i/2rint(^) 
larger than for independent measurements. Stated differently, the number of "ef- 
fectively independent measurements" in a run of length A^ is roughly A^/(2rint(^)). 
The "natural" estimator of ryi(t) is 

1 ^~l*l _ _ 

r^(^) ^ W:rW\ E (^^ - ^)i^^Mt\ - A) . (2.35) 

In order to get a good estimator of rint(^), one sums the terms in eq. (|2.33| ) 
(with r^ computed according to eq. ( |2.35D ) up to |t| < M, where M is a suitably 
chosen cut-off |^. This cut-off is necessary since the "signal" for r^(t)/ryi(0) 
gets lost in the "noise" for |t| ^ r. 

In the following subsection, we describe an alternative method for estimating 
the error, the binning method. Knowing A{A) one can use eq. ( p.34| ) together 
with eq. (|2.3CI|) to estimate the integrated autocorrelation time: 

2.5.2 Binning 

An easy method to analyse the data of a Monte Carlo simulation is the binning 
method. The measurements Ai, i = 1,...,N are first averaged into blocks of 

3 We note that (Ai)mc = (^)- 
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length B called bins 



1 bB 



Ab = J, E ^- b=l,...,NB=[N/B]. (2.37) 

-" j=l+(6-l)B 

If A^ is divisible by B, the average over the blocked measurements is the same as 
the average A in eq. ( |2.28| ). The variance computed from the blocked measure- 
ments is 

_ 1 ^s / 1 ^s \^ 

var(A5) = ttE Ab-t^E^^'.b • (2.38) 

The blocked measurements still suffer from autocorrelations. The error A(y4, i?) 
of the average A, estimated through the variance eq. ( |2.38| ), is 



var(A,5) 



(A(Afi))^ = ^;;/ - (2.39) 

It increases with the bin length B: if the integrated autoccorelation time rint(^) 
is small with respect to B, the systematic effect due to autocorrelations is propor- 
tional to rint(^)/-B |]5S[. The relative statistical uncertainty of the error estimate 



eq. (|2.39|) is approximately given by (2Nb)~^^'^ [^. Increasing the value of B 
the error eq. ( |2.39| ) flattens and oscillates around its correct value ^{A), if the 
number of measurements is large enough to see this. The integrated autocorrela- 



tion time can then be estimated as in eq. (|2.36|) , with Anaivc(^) = A (A, B = 1). 



The error of this estimate is dominated by the uncertainty of A (A) and is given 
by 

A(rint(A)) _ , A(A(Z)) 
rUA) ~ A(A) • ^'-^'^ 

In order to illustrate the binning method, in Fig. |2.2| we show the error esti- 
mate eq. (|2.39| ) as function of the inverse bin length 1/ B for different observables. 



The measurements are performed on a 8^ lattice, for the parameter set (3 = 2.2, 
K, = 0.274 and A = 0.5, after each iteration updating (the Monte Carlo time unit 
is therefore 1 iteration). The observables are the plaquette 

P = ^EE{l-^trP,.(x)}, (2.41) 



the Higgs field length squared 



l^J2^\xMx), (2.42) 
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Figure 2.2: Here, we show the dependence of the error estimates eq. ( |2.39| ) on the 
inverse bin length 1/B. The systematic effects due to autocorrelations, which 
are proportional to Ti^t/B, are clearly visible. The observables are defined in 
eq. (|2.41| ), eq. ( |2.42| ) and eq. ( p.43|) and were measured on a 8"^ lattice with pa- 
rameters P = 2.2, K = 0.274 and A = 0.5. The statistics is 320,000 measurements. 



and the time-likeQ gauge invariant link 
Lt ^ 



-^J22^''[v^i^)Uix,0Mx + 



al 



(2.43) 



We make use of translation invariance and isotropy on the lattice to average the 
observables: this helps to reduce the statistical errors. Reliable error estimates 
for all these observables can be read off at 1/B = 1/400. The number of blocked 
measurements is then Nb = 800 giving a relative statistical uncertainty of the 
error estimates of 2.5%. The systematic uncertainty can be estimated from the 
dotted lines in Fig. |2]^, precisely from the difference between the errors at 1/B = 
1/400 and 1/B = 0, the latter being extrapolated. For all observables represented 
the relative systematic uncertainty of their error estimates is 4%. From eq. (|2.36| ) 
we obtain the integrated autocorrelation times 



Tint(p) 



12.1 



T; 



intl 



10.6, rint(Lt) = 14.2 



(2.44) 



Time and space are the same on a L^ lattice. 



in units of iterations. From eq. ( |2.4CI| ) we can estimate the relative uncertainty 



of the integrated autocorrelation times to be 8%, coming from the systematic 
uncertainty of the errors of the observables. 

The same analysis, repeated for a 32"^ lattice and parameter set (3 = 2.4, k = 
0.2759, A = 0.7, on 24,000 measurements, gives the integrated autocorrelation 
times (in units of iterations) 

rint(p) = 4.9, ri„t(0') = 6.0, T,^tiLt) = 9.9. (2.45) 

The relative uncertainty is again dominated by the systematic effects and is of 
12% for rint(p), 22% for rint(0^) and 14% for rint(-^t)- The autocorrelation times 
have all the same order of magnitude: we conclude that measurements effectuated 
only after 30 iterations should almost be statistically independent. This is in fact 
confirmed by the measurements of the matrix correlation for the static potential 



and the meson mass, see Sect. 4.2. 



2.5.3 Secondary quantities: jackknife binning 

Secondary quantities are defined as 

y = /(A«,A(2),---), (2.46) 

where / is an arbitrary function of the primary quantities A^^\ A^'^\ ■ ■ ■. The 
function / can be complicated, such as the extraction of eigenvalues of a matrix 
correlation, see Sect. p.l|. The best estimate of a secondary quantity is 



y = f{AW,A(^),---). (2.47) 

To estimate the statistical error of y one can in principle use the binning method 
described in Sect. p.5.2| : the quantities yb,B = fi^bBj^b.B^' ' ') ^^^ inserted in 



eq. ( |2.38| ) and eq. ( |2.39| ) at the place of Ah^B- The problem in practice, is often 



that the bins are too small (because of the time costs of the measurements) 
and they fluctuate too much around y. This problem can be overcome with the 
method of jackknife binning. 

For the primary quantities, we consider the bins A^^b, b = 1, ..., Nb and build 
the jackknife averages 



\b = T^ 7L^b',B, b = l,...,NB, (2.48) 



obtained by omitting a single bin in all possible ways. The index B means that 
A^^B is the complement of the bin Ab^B- Evaluating the secondary quantity y 
with the jackknife averages eq. (|2.48|) we obtain the jackknife estimators 



yb,B = /(4i41---), (2.49) 
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with an average 



Vb 



N, 



J2yb. 



(2.50) 



B b 



The error estimate for y can be obtained from |40, 56 



i^iy,B)f 



(iv. 



B 



{^^Jlvh-m^). (2,51) 



For a primary quantity y = A, eq. ( p.51| ) reproduces eq. (|2.39|) . The error estimate 
eq. ( 2.51 ) can be studied under variation of the bin length B as in Sect 



2.5.2. 



Increasing B, the error estimate flattens and oscillates around the correct error. 
The integrated autocorrelation time for the secondary quantity y can then be 
estimated as in eq. ( p.36| ), the naive error being the error eq. ( |2.51| ) for B = 1. 
The jackknife error analysis is our standard method for estimating statistical 



errors. 
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Chapter 3 
Static-light mesons 



As already described in the introduction, we expect the static potential Vq to 
be described in terms of a pair of weakly interacting static-light mesons at large 
separations r of the static charges. A static-light meson is a bound state of a 
static charge and the dynamical Higgs field. The interaction between two such 
mesons is expected to be of Yukawa-type, mediated by the exchange of light color 
singlet bound states of Higgs and gauge fields. We denote by /x the mass of one 
static-light meson: the static potential is expected to reach the asymptotic (in r) 
value 

lim\/oW = 2/i. (3.1) 

In the Hamiltonian formalism explained in Appendix 0, the static-light mesons 
live in the sector of the Hilbert space with one static charge in the fundamental 
representation of the gauge group. We denote by Of^{x) a set of operators labelled 
by i that, when applied to the vacuum state |0), create meson-type states 

\^) = 6f(x)|0) (2 = 1,2,3,...) (3.2) 

localised around the position x of the static charge. These operators carry a colorQ 
index a = 1,2 and transform, under gauge transformation defined in eq. ( [B.8| ) 
and eq. (|B.9|) , as 



RKA)[OrmaR{A) = AL,(x)[6f(f)],,, (3.3) 

where R(A) denote the operator representation of the gauge transformation {A(a;) G 
SU(2)}. The transfer matrix in the sector with one static charge has a spectral 
representation 

TTq = ElVe-^^'-^'KI, (3.4) 



^ Color is the quantum number associated with the gauge group. 
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where |nq, a) is a basis of eigenstates of the Hamilton operator IH in this sector 
with energies E^"* independent of the color a (the sum over the color multiplicity 
of the states is implicit in eq. (|3.4|)). The mass /i of a static-light meson is defined 



as the difference between the ground state energy Eq and the vacuum energy 



f(o) 



M 77^(0) 



/i = E)^>-Er. (3.5) 

It can be extracted from the correlation 

C^{t) = ixr (E[Of (a^)l5Tf [Of (f)],Tr*)/'^) , (3.6) 

where TTq = J2no I'^o) 6xp(— a£',*^°))(no| is the transfer matrix in the zero charge 
sector (see Sect. p.3.2|) , Z = Tr(TrQ'") is the partition function and T the physical 



time extension of the lattice. In the limits 

(T-t)(Ej°)-Ef)>l and t(Ei'') - 4''^) > 1 , (3.7) 

the correlation in eq. ( |3.6D has the asymptotic behavior 

Cjf(t) ~ a*aie-"', (3.8) 

where a^ = [ai]ab = {Oq, a\[Of^{x)]b\0) and the trace over the color indices of the 
states and of the meson operators is implicit in eq. (|3.8|) . In analogy with the 
reconstruction theorem proved in Sect. |B.4| , one can show that eq. ( |3.6D can be 
rewritten in the path integral formalism as the expectation value 

C^it) = {[Of{x + mriU{x,x + m)l[Or{x)]b). (3.9) 

The static charge is represented by a straight time-like Wilson line U{x,x + tOY 
connecting x with x + tO. The meson state \i) is represented by the composite 
field Of^{x) involving Higgs and gauge fields at equal time Xq. To any of such 
fields we can uniquely associate a field operator in the Hilbert space by replacing 
the fundamental fields with the multiplicative field operators defined in eq. ( [B.3| ) 
and eq. ( |B.4|) . The operator associated with Of (x) is precisely Of (x). The only 
restriction in the choice of the fields Of (x) is imposed by the transformation 
property under gauge transformation defined in eq. ( |A.6|) and eq. ( |A.11|) : the 



field Of^{x) must be in the fundamental representation of the color gauge group. 
In addition to the local field $(a;) we can choose for Of^{x) linear combinations 
which take into account contributions from the neighboring Higgs fields (smeared 
fields) and also more general composite fields, with the intent to reproduce the 
wave function of the meson. The physical picture is that of a cloud of dynamical 
Higgs and gauge fields surrounding and bound to the static charge. 
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Actually, eq. ( p.9| ) defines a matrix correlation function. In Sect. ^A 



we 



describe a variational method for extracting from Cf^{t) not only the ground 



state meson mass /i, but also the energy spectrum of the excited states. The 
idea behind the method is that it is possible to find, from the basis of states 



\i) defined in eq. ( |3.2| ), linear combinations approximating the eigenstates of the 
Hamiltonian in the sector with one static charge. The success of the variational 
method is therefore based on the "quality" of the basis of states. This fact makes 



the study of smearing operators important and it is the subject of Sect. p73 . 

Static-light mesons in QCD are a good approximation for B-mesons: the mass 
mb of the h quark is large compared to Aqcd ~ 0.2 GeV and in this sense the 
h quark can be considered a heavy quark. The corrections to the static limit 
rrifo -^ oo are of order Aqcd/''^6 and can be computed in the framework of the 



heavy quark effective field theory, see for example references |59, 60 1. 



3.1 Variational method 

From the matrix correlation function Cij{t) in eq. (|3.6|) (we drop the label M), 
constructed with the basis of states \i) defined in eq. (|3.2| ), it is possible to 
extract the energy spectrum in the charge sector of the Hilbert space with one 
static charge, where the static-light mesons "live". We denote the eigenstates of 
the Hamiltonian in this charge sector by \a). These states have at least a two- fold 
degeneracy: they carry a color index a = 1,2 and their energy is independent of 
the color since the Hamiltonian is gauge invariant. Due to gauge invariance of 
the correlation matrix in eq. ( ^.6| ), the color multiplicity is simply factored out. 
In the following therefore, we drop the color indices of the states and operators. 
Moreover, we restrict our considerations to states with spin 0. This restriction 
is implemented in the way the states are constructed, for example the smearing 
procedures that we employ treat each spacial direction in the same way. The 
eigenvalues of the Hamiltonian are discrete because we are on the lattice, so that, 
in summary, the index a = 0, 1, 2, ... labels the energy levels E^^ which we assume 
are not degenerate. Only the energy differences 

iy„ = Ei'i) - Ef , W^ < W^+, (a = 0,1, 2,...), (3.10) 

have a physical meaning. For the eigenstates we choose the normalisation 

{a\a') = 6aa' ■ (3-11) 

Taking the limit T ^ oo in eq. (|3.6| ), we can write the matrix correlation function 



^ii\ ) as 

Q,{t) = j:{0\6f{xy\a){a\6f{x)\0)e-'^-. (3.12) 
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In practice, the limit T — > oo is reached when T(Ei — E'g ) ^ 1, which means 
that T must be larger than the inverse mass gap in the zero charge sector. This 
is always the case for the situations that we consider, as we discuss in Sect. ^.1| . 
For matrices of the type in eq. ( p.l2| ) a general lemma for the extraction of 



the energies W^ has been proved in |^T|. In this reference, a variational method 
is proposed, which is superior to a straightforward application of the lemma. It 
consists in solving the generalised eigenvalue problem: 

Cij{t)Vaj{t,to) = Xa{t,tQ)Cij{to)Vaj{t,to) , Aq > Aq+i , (3.13) 

where to is fixed and small (in practice we use to = 0). The generalised eigen- 
values \a(t,to) are computed as the eigenvalues of the symmetric matrix C = 
C{to)-^/^C{t)Cito)-^/^ and the vectors 

Va,i = [C(toy^%jVa,j(t,to) with Va,iVa',i = Saa' (3.14) 

are the orthonormal eigenvectors of C. The positivity of the transfer matrix 
ensures that C(t) is positive definite for all t. In [^ it is proven that the energies 
Wa are given by the expressions 

aW^ = ln(A,(t-a,to)/A„(t,to)) + 0(e-*^^") , (3.15) 

where AWa = m.m\Wa — W^l- It is expected that, for a good basis of states. 



the coefficients of the higher exponential corrections in eq. ( |3.15| ) are suppressed 



so that the energies can be read off at moderately large values of t from the 
right-hand side of eq. ( |3.15|) . 

From eq. ( |3.9|) and eq. ( |2.19|) it follows that the matrix Cij{t) is real. Taking 



the complex conjugate of eq. ( |3.12| ), one immediately sees that Cij{t) is symmet- 



ric. In a Monte Carlo simulation these properties are satisfied only in the limit of 
infinite statistics. We make use of the reality property and measure in the sim- 
ulation only the real part of the matrix elements. When we analyse the data we 
symmetrise the matrix by hand. The eigenvalues of C* = C{to)~^^'^C{t)C{to)~^^'^ 
are numerically obtained with the Jacobi method for symmetric matrices ||62| . 
The variational method eq. ( p.l3|) and eq. (|3.15| ) is our standard method for 



extracting the energy spectrum. What we have stated here about this method is 
valid for any charge sector of the Hilbert space. One has to start from a basis \i) 
of states belonging to that charge sector, see Sect. [B.3| . The matrix correlation 
Cij{t) corresponds to matrix elements {j\"ir^\i), n = t/a of powers of the transfer 
matrix operator TT projected into the charge sector. How this works in detail, is 



shown in Sect. |B.4| for the sector with a static charge and a static anti-charge in 
the fundamental representation of the gauge group. The energy spectrum in this 
sector, the static potentials, is the main subject of our work. 
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3.2 One-link integral 

Before describing our choice for the meson-type fields, we would hke to discuss 
a feature of the measurement of the matrix correlation function Cij{t) eq. ( |3.9| ) 
which is independent of that choice. As we see from eq. (|3.8|) , the values of the 
matrix elements Cij{t) fall down exponentially for large t. In order to measure 
these values in a Monte Carlo simulation with statistical significance, also the 
variance of the matrix elements should decrease exponentially^ with t. To achieve 



this, a method called "one-link integral" or "multi-hit" has been proposed in [33], 
which has proven successful. 

The general principle is to replace the observable O, for which one wants to 
decrease the statistical error by another one Oi, with the same expectation value 
but much smaller variance. Such an observable Oj is called improved estimator. 
In the case of the matrix correlation function Cij{t), we observe that it depends 
linearly on the time-like links. When measuring Cij{t), we can substitute the 
time-like links by their expectation values in the fixed configuration of the other 
field variables. These expectation values are called one-link integrals.^ For a 
given time-like link U{x, 0) we write the action eq. ( p.3| ) like 

S = --ti{U{x,0)W^{x,0)} + 

terms independent of U{x, 0) , (3.16) 

W{x,0) = V{x,0) + ^ifi{x)ifi^{x + a6), (3.17) 

where V{x, 0) is the sum of the products of links over the six "staples" around 
the link U{x,0) 

3 

V{x,0) = ^{f/(a;,A;)t/(x + afc,0)f/'^(x + aO,A;) + 



fc=i 



U\x - ak, k)U{x - ak, 0)U{x - ak + aO, k)} . (3.18) 



We denote the part of the action depending on U{x, 0) in eq. ( |3.16D by S{U{x, 0)). 



The expectation value of U{x,0), with all the other field variables kept fixed, is 
given by 



U{x, 0) 



/ dt/(x, 0) f/(x, 0) e-^(^(^'0)) 

(3.19) 



/ df/(x, 0) e-^(f^M)) 
hip) W{x,0) 



Hp) y^det(iy(x,0)) 



^ The alternative is an exponential increase of the number of measurements. 

^ In general, the substitution in an observable of links with their one-link integrals can be 
made under the following restrictions: the observable must depend linearly on the links in 
question and no pair of substituted links can belong to the same plaquette. 
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Figure 3.1: Here, we show the smearing procedure 5*1 for the Higgs field (circles), 
defined in eq. (|3.21| ). The lines represent the link connections. 



where p = (3J(\ei{W{x, 0)) and /„ are the modified Bessel functions of the spec- 
ified integer order. The derivation of eq. ( p.l9|) and the numerical evaluation of 
the ratio of Bessel functions is discussed in Appendix 0. 

An exponential decrease of the variance of an observable with the number of 
links that are substituted by their one-link integrals, is reported for example in 



3^]. The observables considered there are Wilson loops of time extent t and the 
number of integrated links is 2t/a. The variance of Wilson loops computed with 
the one-link integrals decay exponentially with t. 



3.3 Meson- type operators 

We studied different bases of meson-type fields Of (x) by measuring in Monte 
Carlo simulations the matrix correlation function Cij{t) defined in eq. ( |3.9| ) and 
computing from it the energy spectrum of the static-light mesons using the vari- 
ational method described in Sect. pTTI. All composite fields Of^{x), constructed 



with field variables taken at equal time xq and transforming under gauge trans- 
formation defined by eq. ( |A.6| ) and eq. ( [A. Ill ) as 

[Or'\x)l = AL(x)[Of(x)],,, (3.20) 

can be considered. Our aim was to find the best field basis for describing the 
ground state of the static-light mesons. For these studies we simulated the SU(2) 
Higgs model on a 20^ lattice with parameters jS = 2.2, k, = 0.274 and A = 0.5. 
This parameter point is in the confinement "phase" of the model. At the end of 
the section we show the results for the mass of the ground and first excited meson 
state for a simulation at /? = 2.4. The measurement of the matrix correlation is 
improved by the use of the one-link integral method described in Sect. |3.2| . 

We first studied a basis containing the fundamental Higgs field $(x) and 
smeared Higgs fields obtained by iterating the application of a smearing operator 
5*1 to the Higgs field. The smearing operator Si is defined as 

5i$(x) = <l>(x)+ Y. U{x,yMy), (3.21) 

|a; — i/|=a 
a:0=!'0 

where U{x, y) is the link connecting y with x, and is schematically represented 
in Fig. ^TI|. The Higgs field ^(x) is substituted by the sum of itself and of the 
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Figure 3.2: Here, we compare the extraction of the mass /i of a static-hght meson 
using different smearing operators defined in eq. ( p.21| ) and eq. ( p.23|) . 



Higgs fields sitting on the nearest neighbor sites (in the same timeshce) parallel- 
transported to X. Iterating the smearing operator Si we obtain smeared liiggs 
fields 



^t^\x) 



^r ^x) 



(3.22) 



where m = 0,1,2, ... denotes the number of smearing iterations and is called the 
smearing level {m = corresponds to the fundamental Higgs field). We measured 
a matrix correlation function with a basis of smeared Higgs fields corresponding 
to smearing levels 0,1 and 2 of 5*1. The result for the ground state extracted 
according to eq. (|3.15|) is shown in Fig. 3^. We were not able to reach a plateau 
for the ratio ln{Xa{t — a) / Xa{t)) within the range oft considered (up to 8 in lattice 
unit). 

We then investigated a larger basis of meson-type fields, defining in particular 
a smearing operator 5*2 as 



S2<^{x) = V{V^x)+V Y. U{x,y)^{y) 



\x — y\ — y/2a 



V Y. U{x,yMy)} 



(3.23) 



\x — y\^ Vo a 
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Figure 3.3: Here, we show the smearing procedure 5*2 for the Higgs field (circles), 



defined in eq. (|3.23| ). The full lines represent the link connections. 



where P$ = <l>/v$t$ and f/(x, y) represents the average over the shortest link 
connections between y and x. This smearing procedure is schematically repre- 



sented in Fig. |3^. Contributions from Higgs fields sitting on the corners of the 
squares and the cube of side length a around x (lying in the same timeslice as 
x) are taken into account. Through iteration of 5*2 we obtain the smeared Higgs 
fields 



$ 



(m). 



S^ $(x) 



(3.24) 



where m = 0,1, 2,... is the smearing level. We considered the following basis of 
meson-type fields 0^{x), i = 1,2,..., 11: 



Of(x) 
Ofix) 



Ofix) 



Ofix) 



Of{x) 
Of{x) 



Ol{x) 



Of,{x) 



V J2 Uix,yMy) 



\x-y\=a 
^0 = !/0 



V Y. U{x,y)^{y) 

V Y. U{x,yMy) 



$r^)(x), 
1 



5,6,7,8, 



(^{x) X - Y{^\^ - cik)U{x - ak, /c)$(x) + 
" fc=i 

^\x)U{x,k)^{x + ak)} , 

1 
'^(^) ^ T^ Y. {Pki{x) + Pki{x - ak) + 

Pki{x - ak- al) + Pki{x - al)} 
$(x) X ($"l'(x)$(x)) . 



(3.25) 
(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 



(3.31) 
(3.32) 

M/ 



The fields Of^(x), i = 1, ..., 8 have been already described above. The field Og (x) 
is constructed from $(a;) by multiplying it with a "cloud" of gauge invariant links. 
The field 0^{x) is ^(x) multiplied with the sum of the plaquettes around x. 
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Figure 3.4: Here, we compare the extraction of the mass /i of a static-hght me- 
son using all the fields eq. ( 3.25| )-eq. ( |3.32|) and only the smeared Higgs fields 
corresponding to smearing levels 1,2,3,4 of 5*2 eq. ( p.23| ). 



Finally, Oy[{x) is $(a;) multiplied with its length squared. In Fig. |3.4], the result 



for the extraction of the mass of a static-light meson using the fields O 



M, 



X) 



to Fig. 3.2 



11 is shown (triangles). Note the enlarged scale on the y-axis as compared 
We obtain a nice plateau already at moderately large values of t. 



The situation remains practically unchanged (also the statistical errors) if we 
remove from the basis all fields except the smeared fields obtained by iterations 
of the smearing operator 5*2. This means that this smearing procedure contains 
all relevant features for describing the ground state which could be obtained by 
using the larger basis. 

When the generalised eigenvalue problem eq. ( p.l3| ) is solved, the optimal 



linear combination of the basis fields 0^{x) describing the ground state can be 

Therefore, 



M, 



X] 



expressed in terms of the components of the vector vq as J2i Vo,i0, 
we call vq the ground state wave function. Using all the fields eq. (|3.25|) -eq. (|3.32| ) 
for constructing the matrix correlation function, we observe that fo,i, fo,3 and fo,4 
have approximately the same value. This is why we defined 5*2 in eq. (|3.23|) with 
all coefficients in the sum equal to 1. 

Another interesting fact we can learn from the ground state wave function vq, 
is that the field O^, with nearest neighbor contributions, has the lowest coefficient 
vo,2- This explains our original difficulties in extracting the meson ground state. 
In Fig. p.2|, a direct comparison of the smearing operators 5*1 and 5*2, shows clearly 
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Figure 3.5: Here, we show the extraction of the mass of the ground and first ex- 
cited meson state at /3 = 2.4. The basis of meson-type fields was obtained using 
the smearing procedure 5*2 with smearing levels m = 1, 3, 5, 7, 10, 15. The simu- 
lation was performed on a 32^ lattice and the statistics is of 800 measurements. 



that the contributions from the excited states are much more suppressed when 
we use 5*2. 



In Fig. |3]^, we show the results for the static-light meson spectrum that we 
obtained for the parameters (3 = 2.4, k, = 0.2759, A = 0.7 (in the confinement 
"phase") on a 32'* lattice. Details about this simulation will be given in Chapt. ^ 
For the measurement of the matrix correlation function we used a basis with the 
six fields 



<^t\x) 



m 



1,3,5,7,10,15 



(3.33) 



obtained by iterating the smearing procedure 5*2. As we will see in Chapt. H, 
the lattice spacing at /3 = 2.4 is reduced by almost a factor two with respect to 
the lattice spacing at (3 = 2.2. Therefore at (3 = 2.4 smeared fields with high 
smearing levels m are expected to play a more important role than at /3 = 2.2. 
This expectation is confirmed by the simulation. In order to determine with 
confidence the static-light meson masses, we plot in Fig. |3.6| the logarithmic 
ratios on the right-hand side of eq. ( p.l5| ) as functions of the correction terms 
exp{—tAW). This enables us to choose the best time t for reading off the masses 
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Figure 3.6: Here, we show the same data as in Fig. ^3] but plotted against the 
correction term exp{—tAW) in eq. ( p. 151) . We use AW = fi* — ft for the ground 
state and AW = fi** - /i* for the first excited state, see (|3.34| ). 



from the logarithmic ratios and to estimate the systematic errors associated with 
this choice. For the mass of the ground state, we must take the largest value 
t/a = 9. For the mass of the first excited state, we can take t/a = 8. In both 
cases, the systematic error^ are of the same magnitude as the statistical errors. 
However, these errors are small. The results for the meson spectrum are 



a^ 



0.517(2), afi* = 0.88(3), afi* 



1.21(9) 



(3.34) 



We note that the convergence of the right-hand side of eq. (|3.15| ) is not so "crit- 
ical" in the case of the static potentials considered in Chapt. |[ 



"^ The systematic errors for the masses are estimated from the difference between the mass 



read off at the chosen value of t and the crossing point of the dotted hnes in Fig. 3.6 with the 

y-axis (t = cxj). 
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Chapter 4 
String breaking 



We now introduce a method, which - as we will demonstrate in the following 
sections - allows to compute the static potential, Vo(r), at all relevant distances 
in the theory with matter fields. Before explaining the details, we would like to 



mention the basic point, which has first been noted by C Michael [34]. Mathe- 
matically, the method is based on the existence of the transfer matrix [^ and 
the fact that it can be employed also when external static sources are present 



see e.g. fS^j). We have already used this fact in Chapt. y for the computation 



of the static-light meson spectrum. 



As we show in detail in Sect. |B.4| , in the path integral a static source at posi- 
tion X, together with an anti-source at position Xr = x + rk, are represented by 
straight time-like Wilson lines fixed at these space-positions. These Wilson lines 
have to be present in any (matrix) correlation function from which one wants 
to compute the potential energy of these charges. The space-like parts of the 
correlation functions, which are again Wilson lines when one considers standard 
Wilson loops, do not determine which intermediate states appear in the spec- 
tral representation of the correlation functions. They do, however, infiuence the 
weight with which different states contribute. For these space-like parts, we there- 
fore use both Wilson lines which will have large overlap with string-type states 
and Higgs fields with a dominant overlap with meson-type states. Combining 
them in a matrix correlation function, the correct linear combination which gives 
the ground state in the presence of charges can be found systematically by the 



variational method described in Sect. 3.1 



Let us now give precise definitions of the correlation functions, which are 
illustrated in Fig. [4.1| . For small values of r or in the pure gauge theory, the 
static potential can be efficiently computed by means of Wilson loops C^^^^{r,t) 
defined as 

( tr [U{X, Xr) U{Xr, Xr + tO) U\x + tO, Xr + tO) U\x, X + tO)]) , (4.1) 

where Xr = x + rk and U{x, y) denotes the product of gauge links along the 
straight line connecting y with x. For distances significantly larger than the 
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Figure 4.1: Here, the correlation functions used to determine the static potential 
are shown, from left to right: Cww, C*wm and Cmm- The lines represent the 
Wilson lines, the filled circles the Higgs field. 

string breaking distance rt, where the relevant states correspond to weakly in- 
teracting mesons of mass /i, we expect that the potential is close to the value 
lim^^oo Vo(r) = 2/i and can be extracted from the correlation function CMM{r,t) 
defined as 

($"f(x + tO)f/"^(x,X + tO)$(x) <!>^(Xr)U{Xr,Xr + t6)<^(Xr+t6)) . (4.2) 

In order to investigate all (and in particular the intermediate) distances, we 
introduce a (real, see eq. ( p.l9|) ) symmetric matrix correlation function Cij{r,t), 
i,i G {W, M} with C^uir.i) given by 



($^(X + to) U\x, X + to) t/(x, Xr) U{Xr, Xr + tO) <l>(x,. + tO)) 



(4.3) 



In Appendix 0, we construct an Hamiltonian formalism for the SU(2) Higgs 
model with which we can derive the results described above. We summarise 
here the main points of the derivation. The state vectors forming the Hilbert 
space of the theory are represented by wave functionals of the fundamental field 
variables $ and U. The Hilbert space is classified in charged sectors according 
to the gauge transformation property of the state vectors: this transformation 
is related by Gauss' law to the presence of external static charges. The physical 
(gauge invariant) states live in the vacuum sector with no static charges. The 
static potential Vo(r) is defined as the energy of the ground state (normalised to 
the vacuum energy) in the sector with a static quark and a static anti-quarkQ 
separated by a distance r. In Sect. [B.2| , we construct a time evolution operator, 
the transfer matrix in the temporal gauge. We prove that it is strictly positive: 
this allows the definition of the Hamiltonian and ensures the reality of the energy 



^ A static quark (anti-quark) is a charge in the (complex conjugate of the) fundamental 
representation of the gauge group. In the case of SU(2) there is no distinction between quark 
and anti-quark. 
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spectrum. The energy levels in a charged sector can be extracted by evaluating 
powers of the appropriate transfer matrix operatoif] between states belonging to 
this sector. These matrix elements can be shown to correspond to expectation 
values in the path integral formalism. The reconstruction of these expectation 
values from the operator expressions is proved in Sect. |B.4| for the sector with 
a pair of static charges. For charged states generated by the operators given in 
eq. (|B.56| ) applied to the vacuum, the results are precisely the expectation values 
eq. (^), eq. (g^) and eq. 



4.1 Matrix correlation 

The states generated by the operators in eq. ( |B.56|) do not have a space extension. 



In a physical picture we expect the string-type states to be a flux tube |6^, ^, 
68| , |69| , [70| , |7T] of gauge fields binding the static charges. To mimic this situation, 
we introduce smeared gauge fields. For the meson-type states, we expect that 
at large separation of the static charges we can describe the system in terms of 
two weakly interacting mesons. Therefore, we use the one-meson wave functions 
(determined as described in Sect. Of) to construct two- meson states. The one- 



meson states have a space extension due to the smearing of the Higgs field. For a 
high number of smearing iterations, there is effectively an "interaction" between 
the mesons in the two-meson state due to the overlap of the smeared Higgs fields. 
The states entering in the correlation functions for determining the static po- 
tential are restricted by the transformation property under gauge transformation 
and must depend on field variables in the same timeslice. For the string-type 
states we use smeared Wilson lines. They consist of the product of the smeared 
space-like links along the straight line connecting the static charges. We define 
the smearing operator 5* following reference [72 



3 

SU{x,k) = V{U{x,k)+e ^ [U{x,j)U{x + aj,k)U^x + ak,j) + 

\j\x - aj,j)U{x - aj, k)U{x + ak - aj,j)]} , (4.4) 

where V denotes the projection into SU(2). The four space-like "staples" around 
the link U{x, k) are added to it with a weight e which is set to the numerical 
value e = 1/4 and the sum is projected back into SU(2). The smeared space-like 
links corresponding to a number m of smearing iterations are given by 

U^^\x,k) = S"'U{x,k). (4.5) 

For the meson-type states we use the following construction. We determine 



the spectrum of the static-light mesons, using the variational method of Sect. 3.1 



^ The transfer matrix operator in the temporal gauge is restricted to a specific charged 
sector by muhiplying it with the projection operator into the sector. 
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from the matrix correlation function in eq. ( |3.9| ) constructed with the field basis 
Of (x) = $^"'Ha;) [i = 1,2, ..., A^). The smeared Higgs fields $^"'^(a;) are defined 



in eq. (|3.24|) and the numbers rij (i = 1,2, ...,N) denote the smearing levels. 
The eigenvectors Va € H^ (a = 0, 1,2, ...), obtained by solving the generalised 
eigenvalue problem eq. (|3.13|) for large t, are the wave functions describing ap- 
proximately (because of the finite basis of fields and the finite time t) the true 



eigenstates of the Hamiltonian. We define the fields 

N 
i=\ 



^„(x) = Y.^o.,^t\x) (« = 0,1,2,...), (4.6) 



corresponding to the approximate meson eigenstates. The fields we choose as 
basis for the two-meson states are defined as 

\^^{x)\^-\^%{xr)\, a,/3 = 0,1,2, (4.7) 

where x and Xr = x + rk are the positions of the static charges and a,b = 1,2 
are the color indices. The values a = 0, 1, 2 refer to the ground, first and second 



excited one-meson state. The field basis in eq. (f4.7|) contains combinations with 
a ^ (3 which are not symmetric under interchange of the positions x and x^ of the 
static charges. Because we expect the ground two-meson state to be symmetric, 
we project into the symmetric linear combinations of the fields in eq. ( |4.7| ) when 
we analyse the data of the simulations. The "mixed" states (for example of 
one meson in the ground state and one meson in the first excited state) can 
be important when looking at the asymptotic behavior (in r) of excited static 
potentials [^. 

The matrix correlation function, from which the spectrum of the Hamiltonian 
in presence of a pair of static charges can be determined, is constructed with the 
field basis 

mxx)],= [ ^-'(-'"^'^ 1 ^ = 1'2'-'^^ (4 8) 

[UA^.x.)\a. \ [vI^..(x)L [vl/;^(x.)]^ z = iVu + l,...,iVu + 9 ^^-'^ 

where U^™'^\x,Xr) is the product of smeared gauge links (with smearing level 
irii) along the straight line connecting Xr with x and the pairs of indices (ctj = 
0, 1, 2; /5j = 0, 1, 2) label the 9 combinations of meson-type states. Constructing 
correlations like eq. ([4.1|) , eq. ( ^^) and eq. ([4.3|) , but inserting for the space- 
like parts the fields Oi at time xq and Oj at time xq + tO, we obtain a matrix 



correlation Cij{t,r). Its spectral representation is given in eq. (p.6CI| ). We denote 



the energy levels, called static potentials, by V^(r), a = 0, 1, 2, .... In the notation 
of Appendix 0, Va{r) = E^^\r) — Eq . The corresponding eigenstates of the 
Hamiltonian are denoted by |a) = \a,a,b) with color indices a,b. Taking the 
limit of infinite time extension of the lattice T — ;> oo in eq. ( [B.60| ), we obtain 

C^At,r) = 5:(0|6](r)|a)(a|6,(r)|0)e-*^"«. (4.9) 
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fields 


smearing levels m 


S"f/ (seeeq. (|4D) 


7,10,15 


67* (seeeq. (p.23|)) 


1,3,5,7,10,15 



Table 4.1: Here, we list the smearing levels for the gauge and Higgs fields used 



in the simulation with parameters (3 = 2.4, 



K 



0.2759, A = 0.7. 



The operators Oj(r) correspond to the fields eq. ( [4.8| ). The trace over the color in- 
dices of the states and operators is implicit in eq. ([4.9|) . For fixed separation r, we 
extract from C(t,r) the potentials Va{r) using the variational method described 
in Sect. ^A\ 



4.2 Results at /? = 2.4 



Inspired by the investigations in reference [Q, we decided to simulate the SU(2) 
Higgs model in the confinement "phase" near the phase transition line. At fixed 
(3, the mass /i of a static-light meson decreases with increasing k. However, the 
slope (string tension) of the approximately linear piece of the static potential 
for small distances remains constant near the phase transition [^. Thus, string 
breaking is expected to occur at smaller separations of the static charges for larger 
values of k. 

The first results that we obtained [^, ^ were from a simulation at (3 = 
2.2, K = 0.274, A = 0.5 on a 20^ lattice. We observed string breaking at a 
distance Vb/a ~ 5. We decided then to study the system with a better lattice 
resolution at /3 = 2.4. 

The results that we describe in the following are obtained on a 32^ lattice for 
the parameter set 



P 



2.4, 



0.2759, A = 0.7. 



(4.10) 



The field basis is constructed according to eq. ( |4.8|) from smeared gauge {Nu = 3) 
and Higgs fields, whose smearing parameters are summarised in table Table [4.1| . 
The parameters for the simulation were fixed after some trial runs. 

The simulation was performed on a parallel computer CRAY T3E. The 32*^ 
lattice is partitioned in the xy-plane on 8 x 8 processors. We started the sim- 
ulation from thermalised field configurations. The matrix correlation function 
Cij{t,r) {i,j = 1,2, ..., 12) is measured in the ^t-plane starting from each point 
of the lattice up to rmax = 15a and tmax = 9a. The time-like links are replaced by 
their one-link integrals (see Sect. |3.2| ). At r = 1 this replacement is only possible 
for one time-like Wilson line. The matrix correlation function is measured every 
30 iterations of updating (we recall that one iteration updating is composed by 
one heatbath sweep for both gauge and Higgs field, one over-relaxation sweep 
for the gauge field and three over- relaxation sweeps for the Higgs field): this 
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Figure 4.2: Here, we compare the static potential computed from eq. (|4.11 



using the full matrix correlation function (circles) and only the sub-block with 
the (smeared) Wilson loops (triangles). Two representative values of r are shown. 



choice was motivated by the values of the integrated autocorrelation times given 
in (|2.45|) . The CPU time cost of 30 iterations is 80 seconds per processor, for one 
measurement of the matrix correlation function 260 seconds per processor. We 
collected a statistic of 800 measurements. When we analyse the data, we project 
the matrix correlation function Cij{t,r) into the symmetric linear combinations 
of the two-meson fields, thereby reducing its dimension from 12 to 9. 

Autocorrelations in the measurements of the matrix correlation functions for 
the static potentials and the static-light meson spectrum are practically absent: 
the statistical errors, computed by a jackknife analysis (see Sect. |2.5.3|) , remain 
constant when we group the measurements in bins of length 1,2 or 4. 



4.2.1 Static potential Vq 



0, 1, 2, ...) are extracted from the matrix corre- 



The static potentials V^(r) (a 

lation function Cij(t,r) using the variational method described in Sect. pTT| . We 

rewrite eq. (|3.15|) as 



aVJr) 



ln(A,(t - a,to)/Xa{t,to)) + O (e-'^^'^M 



(4.11) 
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where AVa{r) = min |VQ,(r) — Vfl(r)| and the eigenvalues \a(t, to) are obtained by 



solving the generalised eigenvalue problem eq. ( |3.13| ) with the matrix correlation 
function at fixed r. We choose to = 0. 

At all separations of the static charges we compute the static potential Vo(r) 
using the full 9x9 matrix correlation function. As an example, the convergence 
of eq. ( |4.11 ) for r = 4a and r = 11a is shown in Fig. ^]2| (circles). At all distances 



r we can read off with confidence and very good statistical precision (per mille 
level) values for the potential at t = 7a which agree fully with t = 6a. We 
compare these results with what we obtain by considering only the 3x3 sub- 
block of the matrix correlation function corresponding to the (smeared) Wilson 
loops. The resulting potential estimates (triangles in Fig. ^^ are very good at 
short distances but have large correction terms at long distances. Without a very 
careful analysis one might extract a potential which is too high at large distances, 
when one uses the Wilson loops alone. 

4.2.2 Scale ro 

If we want to compute a dimensionful quantity in a lattice gauge theory simula- 
tion, we get a dimensionless number expressing this quantity in units of the lattice 
spacing a as a result. Therefore, we need to fix one dimensionful quantity to its 
physical value in order to get the value of the overall scale a of the simulation. 

In a pure SU(A^) lattice gauge theory, there is only one bare parameter, the 
gauge coupling constant g (or equivalently /3 = 2N/ g"^). In the vicinity of the 
continuum limit, the relation with the scale a is given by the perturbative renor- 
malisation group: 

a = -^ e-^/(26o<;^) (^^^2)-6,/(26g) |^ ^ q(^2) ^ q(^2)| ^ (4^^2) 

llA^ , MN^ 

^°=48^' ^^ = 3(T6^- ('•'') 

The coefficients 6o; &i are the universal one- and two- loop coefficients of the beta 
function. As a result of eq. ( [4.12| ), we see that the continuum limit is reached 
when (7 — > 0. The solution of the renormalisation group equation introduces an 
integration constant Al (called the lattice A-parameter) with the dimension of a 
mass. The development of a dimensionful scale in a theory, which at the classical 
level does not contain any scale, is called dimensional transmutation. 

The eq. ( |4.12| ) is not useful to set the scale a of a lattice gauge theory simu- 
lation because of the 0{g^) corrections on the right-hand side. An efficient and 
precise way of doing it is described in reference fl^ and is based on the force 



between static quarks. A system of two static quarks is approximately realised in 
nature in the cc and bb bound states. The spectra of states of the J/ip and T sys- 
tems are found to be well described by means of a single effective non-relativistic 
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potential |]75[. There are a number of successful potential models (references are 
given in fT^). In lattice QCD we can compute the static potential Vo(r) and from 
it the static force F{r) = dVo(r)/dr. The distance vq is defined through 

r^F{r)\r=ro = 1.65 (4.14) 

and in the phenomenological potentials corresponds to the value 

ro ~ 0.5 fm. (4.15) 

The scale a of a lattice QCD simulation can be set by computing the static force 



and solving eq. ( [4.14| ) to obtain the value of tq in lattice units. Although the 
phenomenological interpretation of this scale is valid only for QCD, the static 
force can be computed in any lattice gauge theory and eq. ( [4.14|) has a solution 
provided the distance rb, at which the gauge string breaks, is larger than tq. Due 
to the clean definition of tq and the good statistical precision with which it can 
be computed, results in lattice gauge theories are often quoted in this unit. 
In order to solve eq. ( |4.14| ) using lattice measurements of the static potential. 



we have first to define the static force on the lattice. We follow ]T^ and define 

aF{rj) = Voir) -Voir -a), (4.16) 

where the argument rj is chosen such that in perturbation theory we have 

nn) = l£^^+Oig^). (4.17) 

To lowest order perturbation theory, the lattice artifacts are exactly eliminated: 
they remain (probably quantitatively reduced) only in the higher Oig^) terms. 
The force defined as in eq. (|4.16|) is called a tree-level improved observable. 

To solve eq. ( [4.14| ) we need to interpolate the force, which is known only for 
discrete values r/a. The general form for our interpolations is 

r'Fir) = f,r~^ + f, + f,r + f,r\ (4.18) 

which corresponds to the potential 

rVir) = -^r~2-/i + /2rlog(r) + /3r^ (4.19) 



The term with coefficient /i is the Coulomb term and the coefficient /s cor- 
responds to the "string tension" (linear term in the potential). To check for 
systematic errors we used three interpolations: (A) two-point interpolation with 
/o = /2 = 0, (B) three-point interpolation with /2 = and (C) three-point inter- 
polation with /o = 0. With the coefficients of the interpolations determined, we 
evaluate the expression for r^/a obtained by solving eq. ( [4.141 ). If this value lies 
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Figure 4.3: Here, the renormalised ground state and first excited state static 
potentials in units of tq are shown as functions of the separation of the static 
charges. String breaking is clearly visible at rb ~ 1.9ro together with the crossing 
of the energy levels. 

in the interval defined by the interpolation points it will be the solution of the 
equation. We want to point out that r^ F{r) in the theory with matter fields is 
not monotonic (because of string breaking we have linij.^oo F{r) = 0): we expect 
that there are two solutions for tq and the smaller one is to be selected. One more 
comment about eq. ( [4.16|) : for Vo(r) we used the values of the ratio Ao(t — a)/Ao(t) 
in eq. (|4.11| ) and repeated the computation of the force and Vq for three values 
t/a = 6, 7, 8. The results for t/a = 8 agree fully with t/a = 7 and are quoted in 
the following. 

For the parameter set in ( |4.10| ) we obtain from all three interpolations (A), 
(B) and (C) the result 



ro/a = 5.29(6). 



(4.20) 



Comparing this number with the values of tq/g computed in quenched QCD [ [76[1 , 
we see that /? = 2.4 in the SU(2) Higgs model corresponds to /3 ~ 6 in QCD. 

4.2.3 Renormalised static potentials Vq [V(r) — 2/i] 

The renormalisation of Wilson loops ly [C] (C is the contour of the loop) in the 
continuum pure gauge theory is considered in references |]7^ |7S[. For a smooth 
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contour C, it is shown that self-energy graphs diverge hnearly in the cut-off A 
with a coefficient proportional to the length L{C) of the contour. The divergence 
originate from the space-time integration region in which the vertices of the graph 
are close together. These divergencies are present in all orders of perturbation 
theory and can be exponentiated 

W[C] = e-'^^^(^) X Wren[C] , (4.21) 

where c is a number ~ 1 and PVren[C*] is a finite function of the renormalised 
gauge coupling. 

In the lattice regularisation, the contributions of self-energies of Wilson lines 
diverge in the continuum like -. We are interested in extracting the static po- 
tentials Va{r) in eq. ( |4.9|) . Therefore, we only have to worry about the diver- 
gent contributions arising from the time-like Wilson lines representing the static 
charges. From the considerations in the continuum, we expect that they expo- 
nentiate with a coefficient proportional to 2t. The same divergencies affect the 
correlation eq. ( p.9|) for the static-light meson and exponentiate with a coefficient 
proportional to t. Therefore, we expect that the quantity 

a [V{r) - 2fi] , (4.22) 

where /x is the (unrenormalised) mass of a static-light meson, is free of divergent 
self-energy contributions and allows the definition of renormalised static poten- 
tials. 

In Fig. ^.3| , we represent the dimensionless potentials vq [V{r) — 2/^] for the 
ground state and the first excited state. For the static potentials aV{r) we take 
the values of the ratios ln(Aa(t — a)/Xait)) (a = 0, 1) in eq. (|4.11|) at large t. The 



computation of the mass /x is discussed at the end of Sect. 3.3. The ground state 



potential shows an approximate linear rise at small distances: around distance 

rb ^ 1.9 ro (4.23) 

the potential flattens. The string breaks! As expected, for large distances the 
potential approaches the asymptotic value 2/i. The first excited potential comes 
very close to the ground state potential around rb and rises linearly at larger 
distances. The scenario of string breaking as a level crossing phenomenon |7S| is 
confirmed beautifully. 

For later purposes, we define a dimensionless renormalised quantity Fi as 

Fi = ro[2/x-K)(ro)]. (4.24) 

The value Vo(ro) was computed using the interpolation eq. ( [4.19|) with three 
parameters /i, /2 and /s (/o = 0). We find the value Fi = 1.26(2). 
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Figure 4.4: Here, we show the overlap ci^o determined from the projected correla- 
tion function eq. ( |4.27| ). The full matrix correlation is used. 



4.2.4 Overlaps 



Overlaps of variationally determined (see Sect. p.l|) wave functions vq are a cer- 
tain measure for the efficiency of a basis of fields used to construct the matrix 
correlation functions. To give a precise definition of the overlap, we define the 
projected correlation function 



n{t) = Vo^iCij{t)Voj = ^Uae 



-tVa,(r) 



(4.25) 



with normalisation! Q{0) = 1 and a labels the states in the sector of the Hilbert 
space with two static charges. The positive coefficients uJa can be derived from 
eq. dU) 



uj^ = |(a|6,(r)|0)|2 with 6,(r) = t;o,.6i(r) 



(4.26) 



and may be interpreted as the overlap of the true eigenstates of the Hamiltonian 
I a) with the approximate ground state characterized by vo- The "overlap" is an 
abbreviation commonly used to denote the ground state overlap, uq. 

We compute vq by solving the generalised eigenvalue problem with C{t = 7a). 
We determine loq straightforwardly from the correlation function Q{t) by noting 



The property r2(io) — vo,iVo,i = 1 follows from cq. (3.14). We use to — 
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Figure 4.5: Here, we show the overlap u^ determined from the Wilson loops 
sub-block of the matrix correlation function. A "naive" way of extracting it 
(triangles), using eq. (|4.27| ), gives an erroneous large overlap at long distances. 
A safe estimate (circles) is obtained from eq. (^.2 



that 



Incuo 



t + a 



In n{t) In n{t + a) (t large) 



(4.27) 



We extract safe values for ujq at t = 7a, which agree fully with t = Qa and are 
shown in Fig. ^^ . Our basis of fields eq. ( |4.8| ) is big (and good) enough such that 
uo exceeds about 60% for all distances. 

It is interesting to consider also the overlap for the Wilson loops alone, i.e. we 
again restrict the matrix correlation function to the 3x3 sub-block associated with 
(smeared) Wilson loops. Let us denote the corresponding projected correlation 
function by fiw(t) and the overlap by u^ . The computation of uj^ is more 



difficult and tricky because it turns out to be very small at large r. In Fig. ^ , 
we present the results for two estimates of u^. The triangles correspond to the 



estimate from eq. (^4.27|) , with Q{t) replaced by i7w(t)- The circles correspond to 
the more reliable estimate using the information from the full matrix correlation: 
the expression 



UJn 



.w 



i-^oo 



Uq- 






(4.28) 
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converges reasonably fast and u^ can be estimated from the r.h.s. for large t. 
Using eq. (^4.28|) , we see that (smeared) Wilson loops alone have an overlap which 
drops at intermediate distances and they are clearly inadequate to extract the 
ground state at large r. On the contrary, using eq. (|4.27|) we get an overlap 
above 50% at large distances: what is estimated here, is actually the coefficient 
uj^ , i.e. the overlap of the (smeared) Wilson loops with the first excited state 
(this statement is supported by direct calculation, see Sect. |4.2.5|) . Because u^ 
turns out to be so large, one should consider Q^^ at much larger values of t 
in order to extract the overlap u^ using eq. ( [4.271) . This might explain the 



problems encountered in QCD for observing string breaking from the analysis of 
a correlation function with Wilson loops only. 

4.2.5 Level crossing 

Finally, we want to get an insight into the interplay between "string states" and 
"two-meson states" in the string breaking phenomenon. The results shown in 
Fig. ^^ support the idea of crossing between the energy levels associated with 
these states. We try to quantify this statement. 

We consider the diagonal sub-blocks of the matrix correlation function eq. ( |4.9| ) 
corresponding to string-type states (fields i = 1,2,3 in eq. (|478|) ) and to meson- 
type states (fields i = 4, 5, ..., 12 in eq. ([4. 8])) separately. We solve the generalised 
eigenvalue problem eq. ( |3.13| ) separately with these restricted matrix correlation 
functions for fixed r and determine approximate ground state wave functions v^ 
for the string-type states and v^ for the meson-type states. With the help of 
these wave functions we construct a projected matrix correlation function 



^ki{t) = <a,(t)^ = ^(V'da)(a|^fc)e-*^'^W {k,l = W, M) , (4.29) 

a 

where 

li^k) = «a(r)) |0) . (4.30) 



In eq. ( |4.30|) , the string-type {k = W) and meson-type {k = M) states are defined 



in terms of the operators that create them when applied to the vacuum |0). The 
definitions in eq. ( [4.30| ) follow directly from eq. ( |4.29| ) and eq. 



Taking t = in eq. (|4.29| ), we get 

nkiio) = iiJim- (4.31) 



The projected matrix correlation in eq. ( [4. 291 ) at time t = is equivalent to the 



scalar product of the states defined in eq. ([4.30|) . From eq. ( p.l4|) , one convinces 



himself that the normalisation of the states {4'k\4'k) = 1 {k = W, M) is a direct 
consequence of the solution of the generalised eigenvalue problem eq. (|3.13| ) with 



to = 0. The scalar product {ipul'ipw) is represented in Fig. ^^ as a function of 
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Figure 4.6: Here, the scalar product of the string-type and meson-type states 
defined in eq. ( [4.30 ) is shown as function of the separation r of the static charges. 



the separation r of the static charges. We see that string-type and meson-type 
states are orthogonal only for large values of r. 
The coefficients 

uju{a) = {a\^k) {k = W,M), (4.32) 

in the expansion eq. ( [4.29| ), express the overlap of the string-type and meson- 
type states with the true eigenstates of the Hamiltonian. The matrix flki{t) is 
real, which means that the coefficients ww(a) and ujm{c() have the same complex 
phase that can be absorbed into a redefinition of the state |a). Therefore, we 
assume that the coefficients uJk{<y) are real. Moreover, we can choose the sign 
conventions co'w(O) > and co'w(l) > 0. We truncate the sum in eq. ( [4. 291) after 
a = 1 and consider the diagonal matrix elements flkk{t) for two fixed times t = ti 
and t = t2'- inserting the known values for Vo(r) and Vi{r), we get a linear system 
of equations for a;^(0) and ujI{1). The solutions read 



^.'(0) 



fifefc(ti)e*^^^W-n,fc(t2)e*^^^W 



ool(l) 



{k- 
{k 



W,M), 

:W,M) 



(4.33) 
(4.34) 



etiAV(r) _ Qt2AV{r) 

e*^^^W[fife,(ti)-^2^0)e-*^^"M] 

where AV{r) = Vi{r) — Vo(r). The sign of the coefficients ti;M(0) and ci;m(1) 
is not fixed yet. From the solutions eq. ( [4. 33] ) and eq. ( |4.34| ), we can compute 
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Figure 4.7: Here, the overlaps of the string-type (circles) and meson-type (trian- 
gles) states, defined in eq. (|4.3CI| ), with the ground state of the Hamiltonian are 
shown as functions of the separation r of the static charges. 



the off-diagonal matrix elements f2wM(^i) and ^wM(t2) using eq. (|4.29|) for the 
four different sign combinations. Comparing with the values that we get from 
the simulation, we can establish the right sign combination. We find that for all 
r, ci;m(0) > and co'm(I) < (in our sign convention). The overlaps eq. (|4.33| ) 
and eq. ( |4.34] ) of the string-type (circles) and meson-type (triangles) states with 
the ground state of the Hamiltonian are shown in Fig. ^^ and with the first 
excited eigenstate of the Hamiltonian in Fig. [4.8| . The results correspond to the 
choice ti/a = 5 and t2/a = 7. Other choices of ti and ^2 give results which are 
compatible within the statistical errors. String-type states have a large overlap 
at short distances with the ground state and at large distances with the first 
excited state. Meson-type states have a large overlap at short distances with the 
first excited state and at large distances with the ground state. In addition, we 
observe that the overlap of the meson-type states with the ground state is also 
large at very short distances. The explanation for this fact is found by looking 
at Fig. ^^, which clearly shows that string-type and meson-type states have an 
overlap with each other at short distances. In the string breaking region around 
r/a = 9 — 10, the overlaps of the string-type and meson-type states have similar 
magnitude, both when the ground state or the first excited state is considered. 
This fact is refiected in the crossing of the energy levels Fig. [4.3|. Here, we would 
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Figure 4.8: Here, the overlaps of the string-type (circles) and meson-type (trian- 
gles) states, defined in eq. ( [4.30| ), with the first excited eigenstate of the Hamil- 
tonian are shown as functions of the separation r of the static charges. 



like to point out that the overlaps represented in Fig. 14.7] and Fig. |4.8| are not 
quantities which have a continuum limit. They are specific to the /3- value and 
the other parameters (e.g. of the smearing) that we consider. 
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Chapter 5 
Scaling 



The lattice spacing a is mainly determined by the choice of the parameter j3 = 
4/(yf^, where g is the gauge coupling. Dimensionless physical quantitie^ F such 
as Fi = tq [2/i — Vo(ro)] contain a dependence on the value of the lattice spacing 
which vanishes in the continuum limit. For a scalar theory we can write [BC, Bl] 

F{a) = F(0) + O ((a/ro)2) , (5.1) 

where O((a/ro)^) summarises terms that contain at least two powers^ of a/ro and 
may be modified by logarithmic corrections. When the corrections in eq. ( |5.1| ) 
(called lattice artifacts) become so small that F is almost independent of a, we 
call this fact scaling. In order to investigate the presence of lattice artifacts in a 
dimensionless physical quantity, it is not sufficient to change the value of /? in a 
model with three bare parameters such as the SU(2) Higgs model. The physics of 
the model is influenced by the choice of all three parameters, for each observable 
differently. For an estimate of the correction terms in eq. ( p.l| ) we must vary 
the lattice spacing a (by changing j3) and tune the bare parameters n and A to 
keep two dimensionless physical quantities Fi and F2 constant. This procedure 
corresponds to the renormalisation of k and A and defines in the parameter space 
a so called Line of Constant Physics (LCP) characterised by 

Fj(/5, K, A) = constant, 2=1,2. (5.2) 



The situation is schematically represented in Fig. |5.1| . Other dimensionless phys- 
ical quantities F3, F4, ■ ■ ■ will, in principle, show correction terms as in eq. (|5.1| ). 
This will give a measure for the scaling behavior of the theory in the investigated 
range of the lattice spacing. 

In Chapt. 1^ and Chapt. ^ we presented results for the static-light meson 
spectrum and the static potentials obtained for the parameter set /3 = 2.4, k, = 



^ We mean ratios of physical quantities with the same mass dimension. 
^ The conjecture that O(a/ro) corrections in eq. (5.1) are absent for any observable in a 
scalar theory, is well accepted. 
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Figure 5.1: Here, a line of constant physics (LCP) in the SU(2) Higgs model is 
shown. The lattice spacing is changed along the LCP keeping two dimensionless 
physical quantities Fi and F2 at fixed values. 



0.2759, A = 0.7 on a 32^ lattice. We want to reproduce the physical situation on 
a coarser lattice corresponding to the choice (3 = 2.2. A sensible choice for the 
dimensionless physical quantities to keep constant, would be to find Fi strongly 
dependent on k, and F2 strongly dependent on A. We have already found the 
right quantity Fi = tq [2/i — Vo(ro)] eq. ( [4.24| ): it mainly depends on the value of 
the mass of the dynamical Higgs field which is in turn determined by the choice 
of K. We now have to face the problem of finding F2. 

Physics shows a dependence on the physical size of the system and on the 
boundary conditions. We regularise the SU(2) Higgs model on a periodic lattice, 
which corresponds to a torus in the continuum and impose periodic boundary 
conditions on the fields. The physical lattice size L/tq is part of the definition 
of a dimensionless physical quantity F. The variation of the lattice spacing in 
eq. ( ^.1|) must be accompanied by a change in the number of lattice points to keep 
the physical size of the torus constant. The dependence of a physical quantity 
on the size of the torus is called finite size effect . 

It is well accepted - supported by the weak gauge coupling expansion p2 
and by early numerical simulations 



m, H - that the SU(2) Higgs model is a 
trivial theory ||3^, which means that the continuum limit is a free field theory. 
Nevertheless, the model can exhibit in a large range of values of the lattice spacing 
scaling properties of a non-trivial almost continuum theory. The interpretation of 
such a behavior is that the SU(2) Higgs model in this range describes an effective 
low-energy field theory and lattice results are perfectly relevant for a continuum 
Higgs model. 
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5.1 Matching of i^ 

Our choice of the first physical quantity to use for the matching of the bare 
parameters along the LCPs is Fi = tq [2/i — Vo(ro)] eq. ( [4.24| ). The first matching 
condition reads 

Fi = ro [2/i - Voiro)] = F* = 1.26 , (5.3) 

where the numerical value is found at /5 = 2.4 with a statistical error AF^ = 0.02. 
Since Fi is mainly sensitive to the value of the meson mass /i, which in turn is 
determined by the mass of the dynamical Higgs field, the quantity Fi is surely 
sensitive to the parameter k. For different values of A, we match the parameter 
K at /? = 2.2 by computing Fi on a 20^ lattice for two values k = ki (giving F^ 
with statistical error AF| ) and k = K2 (giving F| with statistical error AFi ). 
From a linear interpolation 

Fi{ti)\x = ao{X) + K ■ ai{X) , (5.4) 

we determine k* such that Fi{n*) = F^: 

ft:i(F*-Ff) + «:2(Ff)-F*) 



^(1) _ ^(2) 



(5.5) 



Propagating the independent errors AF^ , AF^ and AF*, we get the error for 
the value k* obtained from eq. ( ^.5|) . The results for A=0.5, 0.55, 0.7, 0.757 are 
shown in Fig. p^ As a by-product of these simulations, we obtain the typical 
value of ro/a at /3 = 2.2, for Fi values near Fj*: 

ro/a ^ 2.8. (5.6) 



Using eq. ( [4. 201) , we can derive the change in the lattice spacing between (3 = 2.2 

1.9. (5.7) 



and (3 = 2 A: 

a{p = 2.2) 



a{f3 = 2.4) 

The choice of a 20^ lattice at /5 = 2.2 and a 32^ lattice at /3 = 2.4 corresponds 
only approximately to constant physical volume: we have L/vq ^ 7 at /3 = 2.2 
and L/rQ ^ 6 at /3 = 2.4. But these lattices can be considered relatively large 
for the extraction of the meson spectrum and the static potential: the relevant 



quantity, as we discussed in Sect. pA\ , is the product Lm-n which must be much 
larger than 1. The Higgs mass mn, that we define as the mass gap in the zero 
charge sector of the Hilbert space, can be estimated in our simulations to be 



rQUiii ^ 1 — 1.5, see Fig. |5.5| . Therefore, we do not expect relevant finite size 
effects. 
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Figure 5.2: Here, the matching oi k, a.t (3 = 2.2 on a 20^ lattice using the condition 

eq. ( |5.3|) is shown for different values of A. The matched value k* is obtained from 

eq. Q. 



5.2 Matching of A 

In this section, we propose dimensionless physical quantities and study their A 
dependence. What we are looking for, are gauge invariant quantities with a well- 
defined continuum limit and sensitive to a variation of the parameter A, once the 



parameter k is matched as described in Sect. 5.1 



In the region of small bare couplings, one can use bare lattice perturbation 
theory to determine the LCPs. By using an expansion for small gauge coupling 
one can derive, at the 1-loop level, the following equations ||82| 



dAo(r) 
dr 



1 



r 



dr 



167r2 

1 
167r2 



9 
32^ 



96Ag + ^g^ - 9Ao/ + 



43 



-g^ + 



(5. 



where r = In(amR) ^ is the logarithm of the inverse of a renormalised mass0 in 
lattice units, g^ = 4/(3 and Aq = A/(4fi;2). The dots stand for terms of higher 
order 0(A3,A2^^Ao^^/). 

The change T2—T1 = ln(ai/a2) in eq. ( p.8[ ) can be determined using eq. (|5.7| ). 
If we then use eq. ( p.8| ), with j3i = 2.2 and Ai = 0.5 as initial conditions, to 
compute the change in X a.t j32 = 2.4, we obtainQthe value A2 = 0.8, which is not 

•^ The terms on the right-hand sides of eq. ( |5.8| ) are universal, any renormahsed mass can 
be taken in the definition of r. The non- universal corrections are proportional to a/rg. 
'^ The value of k is taken to keep F\ — F* . 
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far from the value A = 0.7 that we used for the simulations at /? = 2.4. The use 
of eq. (p.8|) is not justified in this case, because the bare couplings are not small 
enough. Nevertheless, we take this crude estimate to start our investigations and 
we consider at /5 = 2.2 the range 0.5 < A < 0.8. 

5.2.1 Cumulants 

We consider the gauge invariant Higgs field S{x) = $"''(a;)$(x) and the construc- 
tion of a renormalised coupling from its connected p-point (or Green) functions. 
There are some subtleties due to the fact that S{x) is a composite field. 

A composite field (or field operator) is a product of fields (field operators) at 
the same space-time point. The renormalisation of composite fields is a compli- 
cated issue and we refer to textbooks, e.g. |85|, |8^, for a detailed discussion. In 
general, a renormalised composite field Ar(x) is expressed in terms of unrenor- 
malised fields B{x) by 

An{x) = Y.ZabB{x), (5.9) 

B 

where Z^b are the renormalisation constants which depends on the bare couplings 
and on the cut-off, e.g. the lattice spacing a in the lattice regularisation. The 
fields B, with which A can mix under renormalisation, have canonical dimension 
equal to or lower than that of A and the same quantum numbers of A. For 
example, the gauge invariant Higgs field S{x) = $''"(x)$(a;) gets renormalised like 

Sn{x) = Z, + ZsS{x). (5.10) 

The canonical dimension of the field S{x) = <l>^(a;)$(x) is two.Q Besides 1 and 
S{x) itself, there is no other gauge invariant scalar field of dimension two or less, 
with which S{x) can mix under renormalisation. 

The connected p-point (or Green) functions of S{x) are defined by the path 
integral expectation values (see later) 

Gcixi,...,Xp) = {S{xi) ■ ■ ■ S{xp))c. (5.11) 

The renormalised connected p-point functions of S{x), which we denote by Gc,r{xi, 
are related to eq. ( |5.11| ) and eq. ( |5.10| ) by 



Gc,r{xi, ...,Xp) — {Sr{xi) ■ ■ ■ Sr{xp))c — ZgGc{xi, ...,Xp) , (5.12) 

The additive renormalisation Zi of the field S{x) eq. ( |5.10| ) cancels in the defini- 
tion of the connected Green functions. The multiplicative renormalisation factor 
Zs in eq. ( p.lO| ) is chosen so that the renormalised connected p-point functions in 

^ On the lattice, the field <i>(x) is made dimensionless by absorbing a factor aj %/k. With 
mass dimension of a field, we always mean the canonical dimension in the continuum. 
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eq. ( |5.12| ) have a well-defined continuum limit, provided all points xi, ...^Xp are 



kept at non-zero distance from one another [q^. The question is, what happens 
in the continuum limit, if some of the points Xi, ..., Xp coincide? This question can 
be addressed in the continuum with the help of the operator-product expansion 



of Wilson ||88[ , according to which 



(5r(x)S'r(2/)), ~ -. (5.13) 

y^x \x — y\ 



The power of the divergence on the right-hand side of eq. (|5.13| ) corresponds to 



the naive dimensional counting. One can easily convince himself that this is true 
for free fields: for interacting fields, logarithmic correction factors may appear. 
We can generalise eq. ( |5.13| ) to the form 



{Sk{.Xi) ■ ■ ■ SB.ixk))c ~ 1 f^ {x ^y). (5.14) 

Xi-^y \x — 2/1^'^ 



We consider now the zero-momentum field 



sr = —jd^xS^^ix), (5.15) 

where V is the physical volume of the torus. The connected p-point functions 
{^r)c {p = 2, 3, 4, ...) can be written as 

(■5r)c = -^ j d^'xi ■ ■ ■ d^Xp (S'r(Xi) ■ ■ ■ Sr{Xp))c 

= :^/dV---dVi('5R(l/i)---M2/p-i)5R(0)),, (5.16) 

where in the second line we used the translation invariance property of the 
connected Green functions. In naive dimensional counting, when yi ^ {i = 
1, ..., p — 1), from the (p — 1) integrations we get 4(j9 — 1) powers of y whereas the 
integrand diverges in this limit like ||/|~^^ {y — > 0), as can be seen from eq. ( |5.14| ). 
Therefore, the p-point functions eq. ( 5.16|) remain finite if 



2p < 4(p-l) ^ p > 2. (5.17) 

Only the 2-point function (s|^)c is naively logarithmic divergent. This problem 
can be cured by defining a modified connected 2-point function 

(4'^ = ^ J dMSn{y)Snm ^sm(^n^ (5.18) 

where T is the time extension of the torus: throughout our work we use T = L = 
V^l^. The sine function avoids contributions coming from y = 0. In order to be 
sure that all divergencies associated with coincident arguments in eq. ( p. 161 ) are 
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regularised for p > 2, it remains to consider the case of two coincident arguments, 
e.g. Hi -^ y2- In contrast with p = 2, there are still two integrations, over yi and 
y2, that regulate the divergent behavior of the integrand, given in this limit by 
eq. ( ^.13| ). We are now able to give a definition of renormalised couplings in the 
continuum 



{S')c 



[{S('^)cf' 



b = 3,4,5,...). 



(5.19) 



The multiplicative renormalisation of the field S{x) cancels in the ratio and we 
can then use the bare fields in eq. (|5.19|) . We call the couplings defined by 
eq. ( |5.19| ) cumulants. We derive in the following expressions for the connected 
p-point functions of the zero- momentum field s. 

In the path integral formalism the p-point functions G{xi, ...,Xp) = {S{xi) ■ ■ ■ 
S{xp)) are constructed |]89[ from the generating functional 



Z[J] = (exp{/d^xJ(x)S(x)}) 

by functional differentiation with respect to the sources J{x): 

5PZ[J] 



Lj^yXi , .. ., Xpj 



J=0 



6J{xi) ■ ■ ■ 6J{xp) 

Combining eq. ( |5.2UD and eq. ( |5.21| ) we can write 

°° 1 /■ 
Z[J] = y^ — / d^xi ... d'^Xp J{xi) ■ ■ ■ J{xp) G{xi, ...,Xp 



(5.20) 



(5.21) 



(5.22) 



The connected p-point functions Gc{xi, ...,Xp) = {S{xi)---S{xp))c are constructed 
from the generating functional 



W[J] = In Z[J] 



in the same manner 



'^cy-^lj •••) -^p) 



5PW[J] 



j=o 



6J{xi) ■ ■ ■ 6J{xp) 
Combining eq. (|5.23|) and eq. (|5.24|) we can write 

^ — / d^xi ... d'^Xp J{xi) ■ ■ ■ J{xp) Gc{xi, ... 



(5.23) 



(5.24) 



W[J] 



1 -^p) 



(5.25) 



p=0 



Expanding the right-hand side of eq. (|5.23| ) in "powers" of J one can derive 
general relations between Gc and G. From these relations, we obtain that the 
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connected p-point functions of s can be expressed for p = 3, 4, 5 by 

{s'), = {s')-3{s){s') + 2{s)', (5.26) 

(/)e = {s'')-4{s'){s)~3{sY + 12{s'){s)'-6{s)\ (5.27) 
is'), = {s')-5{s^){s)-l0{s'){s')+20{s'){sy 

+30{sY{s) - QO{s^){s)^ + 24(s)^ . (5.28) 

Finally, we give a discretised version of the definition of the cumulants that 
we use on the lattice. The zero- momentum gauge invariant Higgs field is defined 
as 

^ = ^T.S{x), (5.29) 

where Q is the number of lattice points. The cumulants are then defined as in 
eq. (lEg) 



<=. = jj^- (P = 3,4,5,...), (5,30) 

where the modified connected 2-point function {s^'^^)c is given by 

{s^'^)c = ^E(^W^(2/))cSin(^^7r), (5.31) 

with {S{x)S{y))c = {S{x)S{y)) — {S (x)) {S (y)) . The connected p-point functions 
{sP)c {p = 3,4,5) are computed according to eq. ( |5.2(j| ), eq. ( |5.27D and eq. ( |5.28| ) 
from the expectation values {s'^) (g = 1, ..., 5). 

Dependence on the volume 

The cumulants Cp for large values of p are very small: in the Monte Carlo sim- 
ulations we were able to obtain a significant signal up to p = 5. In Fig. pT3| , we 
consider the volume dependence of C3 and C4 for the parameter set [3 = 2.4, n = 
0.2759, A = 0.7. We observe that the cumulants are strongly varying functions of 
the lattice size L: they are finite size couplings. The use of the cumulants for the 
matching of the parameters along the LCPs requires therefore a precise matching 
of the physical volume. We decided to compute the cumulants at /5 = 2.2 on a 
6^ lattice: the matching of the physical volume at j3 = 2.4 using eq. (|5.7|) gives a 
lattice size L/a ^ 11.4. 

Dependence on the parameter A 



In Fig. B]Jand Fig. 5.5, we show the A-dependence of the cumulants computed at 



(3 = 2.2 on a 6'^ lattice, with the parameter k matched such that Fi = F* eq. ( [5.31 ). 
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Figure 5.3: Here, we show the volume dependence of the cumulants C3 and C4 
defined in eq. ( |5.3CI| ). 



In the interval of A-values considered and within our statistical precision, the 
cumulants show a weak A-dependence. Especially C4 seems to be sensitive to A. 
However, for the practical purpose of matching the A parameter, this situation is 
not suitable. As can be seen in Fig. O, the cumulant C4 computed at /? = 2.4 
has a strong L-dependence in the region around the matched lattice size L/a ~ 
11.4. The uncertainty coming from an interpolation to this lattice size and the 
statistical errors of C4 at /3 = 2.2 are the limiting factors for the matching of A. 



5.2.2 Higgs and W-boson correlations 

We describe in the following the construction of renormalised quantities with 
a well-defined continuum limit from correlations of gauge invariant Higgs and 
W-boson fields with zero spacial momentum. 

On each timeslice of the lattice, we construct a gauge invariant (composite) 
Higgs field with zero spacial momentum H{xq) as follows: 

H{xo) = (0 ^$t(xo,x)$(xo,x), (5.32) 

where the sum is over the points x in the timeslice with time coordinate Xq. We 
define the correlation function /h(^) as 

^Y.{H{xo + t)H{xo))c, 



hit) 



(5.33) 



where T denotes the time extent of the lattice: throughout our work we use 
T = L. The connected expectation value of the product of two observables Oi 
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Figure 5.4: Here, we show the A-dependence of the cumulants C3 and C4. The 
parameter k is matched such that Fi = F* eq. (p.3|). 
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Figure 5.5: Here, we show the A-dependence of the cumulant C5 and the effective 
Higgs mass Lmff (L/3). The parameter k is matched such that Fi = F^ eq. 



and O2 is defined as usual by 



(Oi02)-(Oi)(02) 



(5.34) 



The correlation function fu{t) is multiplicatively renormalised with renormal- 
isation factor Zg defined in eq. ( p.lO| ). Because the Higgs fields entering the 
connected two-point function in eq. ( ^.33|) are taken at physical time separation 
t, the problems that we encountered in Sect. |5.2.1| , due to the coincidence of 
the arguments of fields in Green functions, are avoided here. From the periodic 
boundary conditions in time direction, we derive the property fniT — t) = /h(^)- 
The correlation function fn{t) is symmetric with respect to the point t = T/2. 
A gauge invariant W-boson field with zero spacial momentum is defined as 



Wrk(Xo) 



L 



'^-itT{TrVk{Xo,x)) 



(5.35) 
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where Vk{x) {k = 1,2,3) is the gauge invariant hnk defined in eq. ( |2.17] ) and 



Tr (r = 1, 2, 3) are the isospin Pauh matrices. We restrict the Lorentz index of 
the W-boson field to the spacial components k. We define the correlation function 

fw{t) as 

/w(t) = ^Y.T.(^rk{Xo + t)Wrk{Xo))c. (5.36) 

^0 r,k 

The real field Wrk{x) = — itr (T,.Vfc(x)) corresponds in the naive continuum limit 
to the field Re tr {Trip'' {x)Dk'^{x)), where Dk = dk + igAkrix)^ is the gauge 
covariant derivative for the gauge field Akr{x) and (p{x) is the Higgs field in 
2x2 matrix notation. The canonical dimension of the field Wrk{x) is therefore 
three and under renormalisation it can mix with fields of equal or lower canonical 
dimension which are gauge invariant and have an isospin index r and a Lorentz 
index k. Because we are in a scalar theory we need a derivative dk to obtain 
a Lorentz index. The only possibility is the field dktT{Trtf'^{x)(f{x)), but it can 
be rewritten like 2Iie tr {Trif^{x)Dk(f{x)). The absence of mixing with other 
fields implies that the correlation function /w(^) is multiplicatively renormalised. 
Again, we have per construction the symmetry property fw(T — t) = fwit). 
From the correlation functions fi{t), z = H, W we construct the quantities 

^^^{t/T) = ln(-^^), z = H,W. (5.37) 

The multiplicative renormalisation of the correlation functions cancels in the 
ratio and keeping the physical size L = T constant the /ij's have a well-defined 
continuum limit. 

Using the transfer matrix formalism it is easy to show that the correlation 
fnit) behaves like 

,T 



/h(^) ~ const X cosh ( ( t) mu ) , (5.38) 



where mn = E[ — Eq . The energies E^^ (a = 0, 1, 2, ...) are the spectrum in 
the zero charge sector of the Hilbert space. We consider mn to be the "Higgs 
mass" in the confinement phase of the model. The relation in eq. (|5.38|) is valid 
in the limits 

t Amu > 1 and (T - t) Amn > 1 , (5.39) 

where Attih = E2 — Ei . We can define an effective Higgs mass m^(t) through 
the relation 

CO* ((I -*).<«)) = ^. (5.40) 
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The effective mass mff (t) converges for large enough t to the Higgs mass ttih- 
We could think to use Lmff (t), determined at a fixed physical time t, as the 



dimensionless physical quantity F2 for the matching of A. In Fig. |5.5| , we show 
the A-dependence of Lm^{L/3) computed at j3 = 2.2 on a 6^ lattice with the 
parameter n matched such that Fi = F^ eq. ( |5.^ ). The effective mass is essentially 
a flat function of A within the statistical precision and it would not be suitable 
for the matching. The same conclusion is valid if we consider an effective mass 
extracted from the correlation function fw{t)- We note that these observations 



are in agreement with the exploratory results of reference [42 



5.3 Scaling 

The results of Sect. |5.2|, concerning the search for a dimensionless physical quan- 



tity F2 which is sensitive to a variation of the parameter A once the parameter 
K is matched using the condition eq. (|5.3| ), can be summarised in the following 
statements: 



The cumulant C4, defined in eq. ( ^.30|) , shows a A-dependence, as can be 



seen from Fig. |5.4|. However, it is a strongly varying function of the lattice 



size, see Fig. |5.3| . The use of C4 for matching the parameter A requires a 
precise matching of the physical lattice size, which can be done only by 
interpolation and is therefore difficult. 

The Higgs and W-boson effective masses are essentially fiat functions of A, 
see Fig. Ol 



In the following we use the assumption, supported by these results, that the 
physics in the confinement "phase" of the SU(2) Higgs model is weakly dependent 
on the parameter A once a dimensionless physical quantity such as Fi = tq [2/i — 
Vo(ro)] is kept fixed. We compare the results for several dimensionless physical 
quantities computed at /5 = 2.2 and at /? = 2.4 for the parameter sets 



and 



/3 = 2.2, 0.5<A<0.8, Fi(/t) = F* = 1.26 (5.41) 



/5 = 2.4, A = 0.7, /t = 0.2759. (5.42) 



This comparison, by a change in the lattice spacing of almost a factor two 
eq. (|5.7|) , gives a measure of the scaling properties of the SU(2) Higgs model 
in the confinement "phase" . 



59 



0.8 



0.6 



=1 0.4 - 



0.2 



T — I I I I I — I I I I I I I I I I I I I 



I I I I 



o 13 = 2.4, L/ro=2.27 

□ i8 = 2.4, L/ro=1.89 

A ,8 = 2.2, X = 0.5 

X ,8 = 2.2, A = 0.7 



m i 



J I I I I I I I I I I I I I I I I I I I T I I I 



1. 




0.1 0.2 0.3 0.4 0.5 

t/T 

Figure 5.6: Here, we show the scahng of the logarithmic ratio /in defined in 
eq. (|5.37| ). The values at P = 2.2 are computed with lattice size L/vq ^ 2.14. 



5.3.1 Cumulants 

Because of the mentioned strong dependence of the cumulants on the physical 
lattice size and the difficulty to match it precisely between /3 = 2.2 and P = 2.4, 
we can only give some indications in form of the following table. 



/3 = 2.2 

c^{L/rQ = 2.12) = 2.39(3) (A = 0.5) 
C3(L/ro = 2.16) = 2.40(3) (A = 0.7) 



C4(L/ro = 2.12) 
C4(L/ro = 2.16) 

C5(L/ro = 2.12) 
C5(L/ro = 2.16) 



1.77(11) (A = 0.5) 

2.08(13) (A = 0.7) 

-9.3(4) (A = 0.5) 

-7.6(5) (A = 0.7) 



/3 = 2A 

C3(L/ro = 2.08) 
Cs{L/ro = 2.27) 

C4(L/ro = 2.08) 
c^{L/ro = 2.27) 

c,(L/ro = 2.08) 



2.21(8) 
2.84(6) 

1.9(3) 
5.1(3) 

-6.4(9) 



The numbers are compatible but no precise conclusion can be made. 



5.3.2 Higgs and W-boson correlations 



In Fig. |5.6| and Fig. |5.7| we show the results for the logarithmic ratios /in and 
/iw defined in eq. (|5.37]). Because we cannot match exactly the physical size of 
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Figure 5.7: Here, we show the scahng of the logarithmic ratio /iw defined in 
eq. (|5.37| ). The values at /3 = 2.2 are computed with lattice size L/tq ^ 2.14. 



the lattice, which at /3 = 2.2 has the value L/vq ^ 2.14, we show for (3 = 2.4 
the results for two different lattice sizes L/tq = 1.89 and L/tq = 1.21 . Moreover, 
for /3 = 2.2 we choose two values A = 0.5 and A = 0.7. The results show rough 
compatibility with scaling, even more pronounced for /i^. 



5.3.3 Static potentials 



In Fig. ^.8| , we compare the results for the renormalised ground state and first ex- 
cited state static potentials that we obtained at /5 = 2.4 (the same as in Fig. [4.3| ) 
and at /3 = 2.2, k = 0.2737, A = 0.5. The results are compatible with scaling 
within minute errors! In Fig. p.9| , we compare the static potentials for two dif- 
ferent values A = 0.5 and A = 0.7 at /3 = 2.2. For A = 0.7 we use k = 0.2928. 
There is no significant difference, confirming the almost independence on A that 
we observed in Sect. ^.2.2| for the Higgs and W-boson mass. These results are a 
strong indication for a continuum-like behavior of the static potentials already at 
the small P values that we use. 
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Figure 5.8: Here, we show the scahng of the renormahsed ground state and first 
excited state static potentials. 
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Figure 5.9: Here, we show the A-independence of the static potentials once the 
parameter k is matched by the condition Fi = F^. 
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Chapter 6 

Conclusions and Outlook 



String breaking, the flattening of tlie static potential at large distances, is clearly 
observed in our lattice simulations of the four- dimensional SU(2) Higgs model 
in the confinement "phase". We are able to determine the ground state and the 
first excited state static potentials with good precision. Our results confirm the 
interpretation of string breaking as a level crossing phenomenon between "string 
states" and "two-meson states". 

In the path integral formalism the static quarks are represented by straight 
time-like Wilson lines, the string-type states are represented by (smeared) space- 
like Wilson lines and the meson-type states by (smeared) Higgs fields. With 
these ingredients, a matrix correlation function can be constructed from which 
the static potentials are extracted using a variational method. We are able to 
describe the level crossing between the string-type and meson-type states in terms 
of the overlaps of these states with the true eigenstates of the Hamiltonian. In 
order to understand the difficulties in observing string breaking in recent lattice 
QCD simulations with dynamical fermions, we studied the extraction of the static 
potential from the string-type states (Wilson loops) alone. The overlap of the 
string-type states with the true ground state at large separations drops under 
10%. The Wilson loops have in this region a large overlap with the first excited 
state. If we would extract the static potential from the Wilson loops alone we 
would see the continuation of the linear rise of the static potential beyond the 
expected asymptotic value, as observed in the QCD simulations |2^, |2^ . 

We addressed the question of the lattice artifacts in our results. The lines 
of constant physics (LCPs) in the parameter space of the SU(2) Higgs model 
are determined by the values of two dimensionless physical quantities Fi and F2: 
the lattice spacing is varied by changing (3 and the bare parameters k and A are 
renormalised in order to keep Fi and F2 constant. The quantity Fi that we choose 
is defined in eq. ( |4.24] ) : it is mainly sensitive to the mass of the dynamical Higgs 



field and hence to k. We studied then the sensitivity to A of different quantities. 
Our aim was to find a renormalised quartic Higgs coupling. We found a candidate, 
the cumulant C4 defined in eq. ( |5.30D , but no precise conclusions could be drawn. 
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Our results support the conjecture that the physics in the confinement "phase" of 
the SU(2) Higgs model is weakly dependent on A, once a physical quantity such 
as Fi in eq. (4.24) is kept fixed. Assuming this, we studied the scaling behavior 



of different quantities between (3 = 2.2 and (3 = 2.4, for a variation of the lattice 
spacing by almost a factor of two. The ground state and first excited state static 
potentials are compatible with scaling within surprisingly small errors: our main 
result is shown in Fig. |5.8| . Compatibility with scaling is also observed for the 
other quantities that we investigated. 

The extension of the method described in |Q and in Sect. § of our work to 



QCD should present "only" problems related to statistical accuracy. There are 
two disadvantages in QCD with respect to our situation: there is no one-link 
integral to reduce the statistical variance of the correlations at large times and 
the computation of quark propagators is so CPU-time consuming that one cannot 
take advantage of translation invariance to reduce statistical errors as we did. On 
the other hand, in QCD the quark fields are integrated out analytically, which 
usually results in correlation functions with relatively small statistical errors. 
As concerns the computation of the quark propagators, the maximal variance 
reduction method of reference [Q is a promising prospective. 

The search for a method of defining LCPs in the confinement "phase" of the 
SU(2) Higgs model is still an open and interesting question: the basic problem 
behind it, is to find a definition of a gauge invariant renormalised quartic Higgs 
coupling. 

There are other important examples in QCD of phenomena involving mixing 
of states, like the computation of the glueball spectrum with dynamical fermions. 
The improvements that were made in the observation of string breaking can 
inspire the research in these other contexts. 
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Appendix A 
Notation conventions 

A.l Lattice notation 

We consider in this work a four-dimensional hyper-cubic lattice in the Euclidean 
space with lattice spacing a 

r = {x/a G N^ I a < x^ < L^, Lq = T} . (A.l) 

We use Greek symbols such as /x, u, to denote all directions 0,1,2,3 and Roman 
symbols such as k to denote the space-like directions 1,2,3. fi is the unit vector 
in /i direction. The physical lattice volume is denoted hj V = TL1L2L3 and the 
number of lattice points by Q = V/a'^. The Fourier transformation of a function 
f{x) on r is defined by: 

f{p) = a^ ^ e-'^-V(x) , (A.2) 

rcer 

where p ■ x = p^x^ = J^fiPf^^f^- Considering periodic boundary conditions in all 
directions: 

fix + L^fi) = fix), (A.3) 

the allowed lattice momenta are restricted to the first Brillouin zone 

B = {p|p^ = — n^, n^ = 0,l,...,L^-l}. (A.4) 

The inverse Fourier transformation of eq. ( |A.2|) is 

/(^) = 4t E e^"- V~(P) • (A.5) 

■p&B 
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A.2 SU(A^) gauge fields 

A SU(A^) gauge field on the lattice is represented by the set of oriented links 
U{x,fi)ab G SU(A^) {a,b = 1,2, ...,N are the color indices) connecting the points 
X + aft with X. A gauge transformation is defined by a field {A(x) G SU(A^) | x G 
r} and by the transformations 

U^{x,^) = A\x)U{x,fi)A{x + afi). (A.6) 

For the pure gauge theory, we use Wilson's action which is based on the plaquette 
P^iv{x) = U{x^jj)U{x + ajl^v)U\x + ai'^iJt)U\x^v) (A. 7) 

and reads 

-^w = /?EE|l-T7RetrP^,(a:)| . (A.8) 

For vanishing lattice spacing and smooth gauge fields, S'w reduces to the classical 
Yang-Mills action S'ym for gauge fields on the continuum 

^w -^ SYM + 0{a^), (A.9) 

if the bare gauge coupling constant g on the lattice satisfies the relation 

2N 
(^ = -T- (A-10) 

A.3 SU(2) Higgs field 

A SU(2) Higgs field on the lattice is a complex scalar field $a(a;) (a = 1,2) 
in the fundamental representation of the gauge group SU(2): under a gauge 
transformation {A(a;) G SU(2) | a; G F} 

$^(x) = A\x)^[x). (A.ll) 

U{x^ \x) is the parallel transporter for the Higgs field, f/(x, /i)<l>(x + a/i) transforms 
under gauge transformation like $(a;). The lattice covariant derivative is defined 
as 

V/,$(a;) = -[?7(x,/i)$(x + a/i)-$(a;)]. (A.12) 

A peculiarity of the gauge group SU(2) is that the complex conjugate repre- 
sentation 2* is equivalent to the fundamental representation 2: 

[ir^Y^Uiir^) = U*, f/GSU(2), (A.13) 
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where r,. (r = 1,2,3) denote the Pauh matrices. This means that the field 
$(x) = 2X2$* (x) transforms hke $(x). Using this property a representation by a 
2x2 matrix field (p{x) can be constructed for the SU(2) Higgs field. We fix the 
notation 

^' [ <^2ix) = M^) + ^Mx) ) ' ^ ' 

where (j)i{x) {i = 1, 2, 3, 4) are the four real components of $(x), and we define 

^i(x) ^i(x) 



^'W^i*;;;) ;;;;;j = v**w*w-w. (a.i5) 

where a{x) G SU(2). Under gauge transformation <f{x) transforms like 

<f^{x) = A\x)^{x) (A.16) 

The scalar product of two SU(2) doublets $ and \E' can be written using the 
corresponding matrices (f and ip constructed as in eq. ( |A.15| ): 



Re$t^ = -tr((^V), (A.17) 

Im$^^ = -tr((ir3)(/?V)- (A.18) 
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Appendix B 

The transfer matrix in the SU(2) 
Higgs model 



B.l Hamiltonian formalism 



Along the lines of reference |6^ we construct a canonical, Hamiltonian formalism 
for fields living on a timeslice. We consider a three-dimensional lattice with 
sites X = a ■ {ni,n2,ns), n^ G N, 1 < A; < 3, I < Uk < N and assume 
periodic boundary conditions. The Hilbert space Ti. is the set of all complex- 
valued, square-integrable wave functions / with arguments U{x, k) G SU(2) and 
$(x) G C^ (in the following we shall denote them simply by f{U, $) or |/)). The 
scalar product is defined by 



im 



D[U]B[<P]f{U,^rgiU,^) 



(B.l^ 



|D[f/]D[0] = l[JdU{x,k)-Y[JdMS)---dMS), (B.2) 

x,k X 

where dU is the Haar measure on SU(2) and (pi, i = 1,2,3,4 are the four real 
components of $, see eq. ( |A.14|) . The field operators ^a{x), Uab{x,k) act as 
multiplicative operators on the wave functions /: 

'Uab{S,k)f]{U,<l>) = U,k{x,k)f{U,<!>), (B.3) 

$„(f) /]([/,$) = $,(f)/(f/,$). (B.4) 

A gauge transformation is defined by a field of SU(2) matrices A(x): 

U^{x,k) = A^{x)U{x,k)A{x + ak), (B.5) 

: At(f)$(f). (B.6) 



$^ff) 



In the Euclidean framework this corresponds to a time independent gauge trans- 
formation. In 7i we define a unitary operator representation R(A) such that 

R(A)/1 ([/,$) = /(f/^,<|.^) ^ /^(t/,$). (B.7) 



69 



Requiring gauge covariance of eq. ( |B.3| ) and eq. ( [B.4| ) one has 

R^ (A) Uab{x,k) R{A) = Al^,{x)Ua'b'{x,k)Ab'b{x + ak) 
Rt(A)|.,(f)R(A) = AUx)^a'{x). 



An operator O is called gauge invariant if 

R(A), 6 



0. 



(B.8) 
(B.9) 



:b.io) 



In the following we shall need operators O which can be represented with the 
help of an integral kernel Ko{U, $; U', $'): 



(6/)(f/,$) 



D[U'] D[0'] KoiU, $; U', $') /(f/', $') . (B.ll 



The gauge invariance condition eq. ( |B.10|) translates into 

For the product of two operators Oi O2 we have 

KoMU,<^;U',<^') = Jd[U"]D[c^"]KoAU,<^;U",<^" 

The trace^ of an operator is defined as 

TrO = /D[f/]D[</)]iro(f/, $;[/,$). 



(B.12) 



:B.13) 



(B.14) 



B.2 The transfer matrix in the temporal gauge 



We define an operator TTtemp through the following kernel: 

i\:T°'P(f/, $;[/',$') = TH(f/,$)-TG(f/)-5(f/, $;[/',$') 

TG{U')-Tn{U\^'), 

where 



(B.15) 



1 ii PjAx) 

2 ■'^ , 



y '^ X l<i<j<3 

TH(?7,$) = exp|-^5] A($t(f)$(f)-1 



:B.16) 



(B.17) 



^ We write Tr for the operator trace in 7i and tr for the color trace in the gauge group. 
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-kJ2 {^\x)U{x, k)^{x + ak) + ^\x + ak)U\x,k)^{ 



X 



fc=i 



fi;(<l>t(f)$'(f) + $'t(f)$(f; 



^(t/,$;t/',$')=exp|5: 

-- (<^\x)<^{x) + $'"f(x)$'(f)) - /3 ^ ('l - i tr (f/(f , A;)t/'"f(f, A;)) 



(B.18) 



We included the terms quadratic in <l> and $' separately, in S eq. ( |B.18|) instead 
of Th eq. ( |B.17|) for later convenience. 



B.2.1 Posit ivity 

Proposition 1. 

a) TTtemp is a selfadjoint, bounded and gauge invariant operator in Ti. 

b) TTtemp is strictly positive for k > and A > 

(if A = strict positivity holds for < n < 1/6). 

Selfadjointness of i^T""^ means K^^'^iU', $'; U, $)* = K'^^'^''{U, $; U' , $'). This 
follows from the fact that the kernel /f^™'^ is symmetric under interchange of 
{U, $} with {[/', $'} and that the substitutions Uabix) -^ U*i,{x), $a(^) — > "^al^) 
leave the integration measure and the action invariant. 

In order for k!^^"^ to be bounded, one has to discuss the cases A > and A = 0. 
If A > the kernel is bounded for all values of k because of the term quartic in 
the Higgs field. If A = one has to show that the matrix 



B. 



xa,yb 



hAh -^T. [Uab{S, k)6.^^^j^ + Uliy, k)6.^y._^^f^] (B.19) 



k=l 



is strictly positive. This is the same as for fermionic fields |^5[: the restriction 

\n\ < 1/6 has to be imposed. 

Gauge invariance of Kr^^"^ (eq. (|B.12| )) is obvious. 

To prove strict positivity of TTtcmp one has to verify that 

D[U] D[0] D[U'] D[0'] f{U, $)* irr"(f/, $; U' , $') f{U', $') 

= (/iTTten^pl/) > for all / ^ in 7i. (B.20) 

If Kr^^"^ is bounded it is enough to show strict positivity for the kernel S{U, $; U', $') 
because the kernels Tq and Th can be absorbed in a new wave function g = 
TgTh f G H. Furthermore, one has 



S{U, $; U', $') 



n Snmx); $'(f)) ■ n S^Uix, k); U'{x, k)) 
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1 



$"^$' + $'"''$ ($"''$ + $'"''$') 



i--tr(f/f/'t; 



S'h($; $') = exp <5 K 

Sg{U;U') = exp|-/3 

and we are left to prove that 

fdUdU'd^(j)d'^(j)'h{U,<^ySni<^;<^')SG{U;U')h{U',<^') > (B.21) 

for all square integrable nonvanishing functions h{U, $) depending on one link 
variable U and one field variable $. Changing to the matrix notation for the 
Higgs field ip = pa, p >0, aG SU(2) one has 



$t<|,' I ^'1^ 



tr (v^V 



t,.M 



d(/)i(x) ■ ■ ■ d(j)Jx) oc / P dp da 

R4 ' ' Jo Wsu(2) 



(B.22) 
(B.23) 



Expanding 



expJKtr {(p^(p')} 
tr fa'l'a;' 



E 

n=0 



(/.pp')" 



n! 



tr (a'' a] 



we recognise that 

SU(2) which is composed o 

end up with 

'S'h(p, a;p', a') = 



is the trace of the tensor product representation of 
n quark^[|^. Reducing out the tensor product we 



g(^^^^(^^^(.)(^t«,^.^(pW^)/2 



n=0 



n\ 



;B.24) 



z/ labels the set of all inequivalent unitary irreducible representations of SU(2) and 
X^'^^ is the character of the representation u. The coefficient b,y{n) is the number 
of times the representation u occurs when reducing out the tensor product. Since 
all irreducible representations can be obtained in this way, for every u exists a 
value of n for which bi,{n) ^ 0. Similarly, the character expansion for Sq is 



Sg{U-U') 

Defining 

hn{U,a) 
and using the relation 



e-^-Y.CuX^'KUU'^), c,>0 



("OO 

/ dp p^-^"" h{U, pa) e- 
Jo 



■pV2 



dUx^''\VU)x^'''\U^W) = S,,>^x^''\VW) 



(B.25) 



(B.26) 



;B.27) 



^ A quark is a vector transforming according to the fundamental representation of the gauge 
group, an anti-quark to its complex conjugate. In the case of SU(2) they are equivalent. 
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where d^ = x^'^K'^) is the dimension of the representation u, the integral in 



eq. (|B.21|) becomes 

^^C(z/,z/',n) fdVdW fdUdahl{U,a)x^''\UV)x'^'''\a^W) 



vm' n 



(B.28) 



where C(z/, z/',n) = dpd^iCui^^hyi{n)e~^ . We are finally left with a sum of strictly 
positive terms provided k > 0. 

One can show that the condition /t > is also necessary. It is possible to 
construct a function /i(f/, pa) = h{p) such that hn eq. (|B.26|) vanishes for all even 
n and is different from zero for at least one odd no. If k is negative, for this 



function the whole sum eq. ( |B.28| ) is then negative. 



B.3 Charge sectors 

B.3.1 Minkowski space continuum 

The classical Lagrange density for the SU(2) Higgs model in Minkowski space is 

1 



C 

D'^$(a;) 
y($(x)) 



d'^A^ix) - d'-A^.ix) - geatcA^,{x)A:{x) 



(B.29) 



5" + wA^ai^: 



Ta 



$(x), 



m, 



l^^{x)^x) + Xo[^\x)^{x)f 



where the indices a, b, c label the three generators of SU(2) y, r^ being the Pauli 
matrices, and g, ttlq, and Aq are the bare couplings. The real vector fields A^ 
carry a Lorentz index p and a group index a. $ is the complex Higgs SU(2) 
doublet field. The canonical conjugate momenta to the fields A^ and $, $^ are 



n:^(^) 



dC 



TT 



d{doAaf,{x)) 

dC 
~ d{doHx)) 

dC 
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-Fr(x) = 

[D°$(x)]t, 
D°$fa;). 



-Ei p = i-- 
/i = 



1 9 S 
' '"^ (B.30) 



(B.31) 
(B.32) 



E\ is the non-Abelian electric field strength. The generalisation of Gauss' law 



follows from the equations of motion in the Lagrange formalism |91|: 



f 

Ja 



(B.33) 
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j^ is the matter-field current and in particular 

r„ , , + T, 



J a 



-w 



' a. ^ ;^\ 'a- t 
Vr — $ - $T — TT^ 

2 2 



(B.34) 



Performing the Legendre transformation to the Hamilton density one has to 
face the problem that the canonical conjugate momentum to A^^ is zero. One 
can choose a gauge in which the fields A^ vanish identically, the temporal gauge 
^. There is still some freedom left in the choice of the gauge, namely the time 



independent gauge transformations. The canonical quantisation of the theory 
proceed considering wave functionals 



nK{x),^{x)] 



(B.35) 



where the fields are taken at fixed time. In this representation of the states the 
field operators A\{x) and $(a;) act as multiplicative operators and the conjugate 
momenta are represented by functional derivatives 

5 



El(x) 



vr X 



7r^(a;) 



SAai{x) 

. 6 



5^{x) 
. 6 
''6^Hx) 



Under a time independent gauge transformation A{x) 
finitesimal parameters uja{x) the fields transform like: 



(B.36) 
(B.37) 
rB.38) 



1 — iuOa{x)^ with in- 



5k^{x) 



-d'ujaix) + eabc^b{x)Al{x) 



-iuJa{x)^^{x) 



lU, 



a{m\ 



X) 



Ja 



(B.39) 

(B.40) 
(B.41) 



Now, we would like to find the generators G^ of the gauge transformations on 
the wave functionals \1/: 



5k^ 



-- I d^XUJa{x)Ga{x)^ 

9 



(B.42) 



To this aim we write: 
(5a* = / d^x 



(5\I^ 5\& + 6"^ 

6aAI^{x) + 7:i:7^'5A$(f) + Sa^\x)- 



uAai[X 



5^{x) 



6^^{x) 



— J d'xuja{x)\^ [diE'M + geabcAi{x)Ei{x)] * + 
ig [(7r(f)*)-$(x) - $t(f)^(7rt(f)*) 



(B.43) 
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In the second line we used integration by part. From this we conclude that 

G„(f) = D,i:(f)-j°(f), (B.44) 

where the operators on the right hand side are defined by 

D,E:(f) = d,El{x) + geabA{x)El{x) , 

j°(f) = -z^[|.^(x)(|)^^7r,(x)-<|t(f)|^t(^)" . 

We can classify the wave functionals \1/ according to the action of Ga- Gauss' law 

Ga^ = (B.45) 

is equivalent to the gauge invariance of the state \1/. On states transforming like 

Sa^ = -tuJa{x)T^''^<i/ , (B.46) 

where Tj^^^ are the generators of the irreducible representation u of SU(2), the 
operator G acts as 

Gaix) = gT^'^S^'^\x-x'). (B.47) 

Thus wave functionals transforming under gauge transformation according to the 
representation u of A(a;) for some fixed spacial point x express the presence of a 
static external charge at position x in the representation z/. 

B.3.2 Lattice formulation 

Strict positivity of TTtemp allows the definition of an Hamiltonian IH: 

H = --InTtemp (B.48) 

(Jj 

From the properties of TTtemp it follows that IH is selfadjoint (has only real eigen- 
values), bounded from below and gauge invariant. 

The gauge symmetry of IH allows to choose its eigenstates with a well defined 
transformation property under the gauge group. The Hilbert space Ti. can then 
be classified according to the irreducible representations of SU(2) on each point 
of the lattice |P2|. As we saw inspecting Gauss' law in the continuum these 



representations define the charge sectors of 7i. The transfer matrix is restricted 
to a specific charge sector by multiplying it with the projector P onto this sector. 
We give two examples of charge sectors which are of interest for us. 

The zero charge sector (vacuum sector) corresponds to the gauge invariant 
wave functions. We denote by |no) the eigenstates of IH with eigenvalues £'^°'' 
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which form a basis of this sector. The corresponding projection operator IPq = 

Eno \no){no\ is 

(Po/)(f/,$) = jY['^Wiz)f'^iU,^). (B.49) 

z 

The transfer matrix operator in the zero charge sector is TTq = TTtempIPo and has 
the kernel 

ii'^(t/, $;[/',$') = Jl[dW{z)K'^'''''{U,<!>-U'^\<!>'^') (B.50) 

z 

A wave function / in the sector of 7i with a static quark at position x and a 
static anti-quark at position y transforms as 

fab = E ^aa'ix)f,,yAUy) (B.51) 

a',b'=l 

We choose a basis |nqq, a, b) with eigenvalues -E'^'^^-' of IH independent of a, 6. We 
drop in the following the color indices a, b of the states when an implicit sum 
over them is meant. The projection operator '\P{x,y) = J2n - l"-qq)('^qql o^ito this 
sector is 

(lP(f,y)/)(f/,$) = f'[[dW{z)AtTW^{x)tTW{y)f^{U,^). (B.52) 

z 

The transfer matrix operator in this sector is Trq(^)q(j^) = TrtempIP(a;, y) and has 
the kernel 

Kf{U,^;U',^') = fYldW{z)AtTW\x)tTW{y) ■ (B.53) 

z 

ir^''°P(f/,<l>;f/'^',<l>'^'). (B.54) 

B.4 Reconstruction of the Euclidean expecta- 
tion values 

We consider operators 0{x,y)ab which transform under gauge transformation 
according to 

Rt(A)6(x,f)a,R(A) = Ai,,{x)6{x,y)a'b'Ab'b{y). (B.55) 

These operators create out of the vacuum states belonging to the sector with a 
static pair of quarks. Examples are: 

^a{x)^{y) and U{x,y)ab (B.56) 

U{x, y) denotes the product of links along a path connecting y with x. Because 
of eq. ( p.3|) and eq. ( p.4| ) we can directly associate the observables 0{x,y)ab = 
^a{x)^t{y) and U{x,y)ab- 
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B.4.1 Reconstruction theorem 
Proposition 2. 

a) 

Tr(Tf ) = n/df/(x,/i) ■ n/#i(^) ■ ■ ■ #4(0;) exp(-5) = Z 



x,ii 



is the partition function on a four- dimensional lattice with time coordinate 
Xo/a = 1, 2, ..., M and periodic boundary conditions. 



b) 



Tr EB(5^,^lL<J),,(,^A(f,y),X'~*/M /Tr(Tf) 



= ( tr (A(x, y)U{y, y + t6)B{x + t6,y + t6)^t/(x, x + tO)^) ) = Cab 

where xq = yo and A, B are operators of the type of eq. (p.55| ). This is the 
amplitude for the transition from the state A to the state B over a time 
interval t. 

We first inspect the kernel of TTq in eq. (|B.50| ). Identifying W{z) with the 
time-like links f/(z, 0), z = {zo,z), renaming $'(z) — > $(2 + aO) and U'{z,k) — >• 
U{z + aO, k) we see that eq. (|B.18|) , excluding the terms quadratic in $ and $' 



separately, becomes 

SiU,<^;U"^\^"^') = exp{-ASizo,zo + a)}. (B.57) 

AS{zq, zo + a) are the terms in the action which couples variables on the neigh- 
boring timeslices t = zq and t = zq + a. From this we see that the integration 
over the time-like links in the path integral is equivalent to the projection onto 
the gauge invariant sector of Iri. eq. ( |B.57D is also the reason why we called TTtemp 



the transfer matrix in the temporal gauge. The other factors in Kr^"^^ contain 
the pieces of the action depending only on the fields in the timeslices and the 
trace is equivalent to the periodic boundary conditions in the time direction. 
The numerator of the left hand side of b) can be written as 

n / dt/(^, /i) ■ n / dV(^) tr iA{x, y)B{x + t6,y + tof 

2*/" tr U\x + {t- a)0, 0) tr U\x + (t - 2a)0, 0) ■ ■ ■ tr U\x, 0) ■ 

2*/" tr U{y, 0) tr U{y + aO, 0) ■ ■ ■ tr U{y + (t - a)0, 0) ■ 

exp{-S) . (B.58) 
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Because the action and the measure are gauge invariant, for any observable O 
one has {O'^) = (O), where O^ is the gauge transformed observable. Gauge 
transforming eq. ( [B.58|) , integrating over JJ^^^^ J JJ^ d A{wo, tv) and applying 
eq. (p.27| ) one sees that the product of traces in eq. (|B.58| ) is glued into one 
trace, the right hand side of b). 

B.4.2 The static potential 

Using the spectral decomposition To = J2no l^o) exp(— ai?^°))(?7,o| the partition 
function can be written as 

Z = ^exp{-MaEi°)} . (B.59) 

no 

Inserting 'T^q{x),q{y) = Sn - I'^qq) 6xp(— a-E^''^-')(nqq| in the quantum mechanical 
expression for the correlation Cab we have: 

Cab = J Yl (E(^o|Bl>qq)Kq|A,e|no) ) • (B.60) 

no,nqq y b,c / 

exp {-t(i?(^^) - E^^^)} ■ exp {-MaE^^^} 

In the limit of Ma—t ^ l/m[ , where rrii is the mass gap in the gauge invariant 
sector, and t ^ {E[ — Eq )~'^ we have the asymptotic behavior 

Cab ~ ftB «A exp {-t(4"^^ - 4°^)} , (B.61) 

where a^ = (Oqq|A|0) and the trace over the color indices of the states and 
operators is implicit. The quantity 

V = Ej^'^^-Ef (B.62) 

is called the static potential. 



Appendix C 
Updating algorithms 



In this appendix, we give a detailed description of the algorithms that we used 
for the Monte Carlo simulation of the SU(2) Higgs model on the lattice. Before 
describing the algorithms, we want to discuss the generation of random numbers 
which are an essential part of the updating. 

C.l Random Numbers 

Because the updating process is stochastic, one needs to generate a large num- 
ber of independent random numbers. We use the high-quality random number 



generator of M. Llischer [£^, ^, derived from an algorithm originally proposed 
by Marsaglia and Zaman |9^. It generates random floating point numbers dis- 
tributed uniformly in the range [0, 1). From these flat random numbers we have to 
construct random numbers with probability distributions needed in the updating 
of the gauge and Higgs field variables. 

To be clear with the terminology, one speaks in general of random numbers r] 
distributed in an interval [a,b] C R with probability density ^(77). The normali- 
sation of the density is such that /^ p{t]) d?7 = 1. A uniform (or fiat) distribution 
of a random number r] in the interval [a, b] corresponds to p{r]) = 1/(6 — a). The 
probability distribution P{rj) is related to the probability density by 

'JM . ,(„). (c.i) 

If one considers a transformation of variable fj = b{ri) the probability density p{r]) 
of the transformed variable is related to p{ri) by the requirement 

dP{r]) = dP{fj). (C.2) 

First of all, we consider the generation of random numbers fj distributed 
according to the Gaussian density 

Piv) = -^e-^'- (C.3) 
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It is better to consider the generation of a couple {rji, 772) of independent Gaussian 
random numbers. Independent means that 

dP{f]i,f]2) = Pivi)dfiip{f]2)df]2. (C.4) 

With the variable transformations fji = pcos6 and 772 = psin^, where p G [0, 00) 
and 9 G [0, 27r) we can rewrite eq. ( |C4| ) using eq. ( |C.2| ) as 

dP{p,9) = dP{fii,fi2) = -pe-p'dpde 

= d(l- <=-'>(!?)• (C.5) 

This is equivalent to the generation of a couple {i]i,ri2) of independent flat random 
numbers in the range [0, 1) together with the transformations rji = 1 — exp(— p^) 
and 772 = ^/(2vr). The final result is then 



rji = J — \n{l — rji) cos{2'Kri2) and 



?72 = -\/-ln(l -?7i)sin(27r?72). (C.6) 

Gaussian random numbers are needed in the updating of the Higgs field variables, 
see Sect. |(J.3.1| . 

Secondly, we consider the generation of random numbers distributed in [0, cxd) 
according to the density 

2 

P{y) = ^v^exp(-y)e(2/) , (C.7) 



where 

e(x-a;o) = I I f' "^^"^^ (C.8) 

^ ^' [ for X < xo ^ ^ 

is the Heaviside function. This can be done by generating a number a G [0, 00) 
according to the Gaussian density 

p(a) = -^e-^'eia) (C.9) 
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and independently a number b G [0, 00) with the exponential density 

p{b) = e-''e{b). (C.IO) 



To generate a one uses eq. (|C.6|) replacing 2Trri2 with 71172/2 as the argument for 
the trigonometric functions. To generate b one simply takes a fiat distributed 
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random number i] in the range [0, 1) and sets b = — ln(l — 7]). We define the 
variables 

y = a and y = a^ + b (C.ll) 

with ranges y G [0, oo) and y G [0, ^/y). From eq. ( |C.2| ) we get the relation 

dP{y,y) = dP{a,b) = -^e-yQ{y)Q{^ - y)dydy , (C.12) 



which means that the variables y and y are distributed with probability density 
piSiU) = 2/-\/vrexp(— y)0(|/)6(y^ — y). The probability density for y alone is 
obtained by integrating p{y, y) over all possible values of y: 

p{y) = f"^ ^e-yQiy)dy = ^^e-ye{y) . (C.13) 

Jo x/iT Jti 



The random numbers y from eq. (|C.11| ) are needed in the updating of the gauge 
link variables, see Sect. |C.2.1 . 



Now, we are ready to describe the different parts of the HOR algorithm in- 
troduced in Sect. Ol. 



C.2 Updating of the gauge field 

To discuss the updating of the link variable f/(x, /i), we split the action eq. 
into 

S = -^ti{Uix,fx)W^ix,fx)} + 

terms independent of U{x,^) , (C.14) 

2k 
W{x,ij) = V{x,ij) + —(p{x)(p'^{x + aft), (C.15) 

where V{x,fi) is the sum of the products of links over the six "staples" around 
the link U{x, fi) 

V{x,ix) = y^{?7(x, v)U{x + ai>, n)U^{x + afi, v) + 

U\x — ai>, v)U{x — ai>, iJ,)U{x — ai> + aft, u)} . (C.16) 
We denote the part of the action depending on U{x, /i) in eq. ( p.l5| ) by S{U{x, /i)). 



C.2.1 Heatbath 

The heatbath algorithm for the gauge field in the SU(2) Higgs model is a simple 
modification of the algorithm for the pure SU(2) gauge theory |^, |93, ^. We 



drop the arguments (x, /i) of the link to be updated. The new link U' is chosen 
independently of U according to the Boltzmann distribution 

dP(t/') ~ e^t^(^'^'M?7'. (C.17) 



The matrix W defined in eq. ( |(J.15| ) can be written as W^ = Vdet W W, where 



W G SU(2). Using the invariance of the Haar measure we obtain from eq. ( |C.17| ) 

dP{U'W) ~ e5^*^(^')df/', (C.l^ 



where p = /3\/det W. Writing U' in the quaternionic representation 

U' = ao + iajTj , a^ G IR , a^ = a^a^ = 1 , (C.19) 

where Tj {j = 1, 2, 3) are the Pauli matrices, the Haar measure for SU(2) takes 
the form || dU' = ^^(a^ - l)d^a and eq. (lOTsD becomes 



dP{U'W) ~ ^6{a^ - l)eP''° daod^a . (C.20) 



7r2 



The heatbath algorithm consists then in the following: 
1. Generation of a^ according to the distribution 

dP(a^) ~ Jl - ay^daod^nSin^ - 1) , (C.21) 



where aj = rijJl — Oq. This is done using the method described in PB 



which is a slight variation of that of Fabricius and Haan [Q. In order to 
generate Oq G [—1,1] with the distribution 



dP(ao) ~ ^Jl-ay°dao (C.22) 

we perform the change of variable 

y = p(l-ao)G[0,2p]. (C.23) 

According to eq. ( |C.2|) we rewrite eq. ( |C.22|) as 



dP{y) ~ (2-^)V2y^e-M2/. (C.24) 

The generation of random numbers y according to eq. ( |C.24| ) cannot be 
done exactly. Instead, one generates y according to the probability density 
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eq. ( |C.7| ) as explained in Sect. |C.1| . An accept-reject step takes into account 



the omitted factor (2 — ^)^/2. A flat random number r] E [0, 1) is generated: 
if 

27^2 < 2-^ (C.25) 

P 

the change oq = 1 — y/p is accepted. We notice that the values y > 2p are 
automatically rejected in eq. (|C.25|) . The components Uj of the vector n 
in eq. ( |C21D are uniformly distributed on the surface of the unit sphere in 
three dimensions. The generation of rij requires two flat distributed random 
numbers rji and ri2 in the interval [0,1): 

ni = l-2rii, (C.26) 



ns = Jl-nlcos{27rr]2), (C.27) 



ng = Jl-nism{2TTr]2). (C.28) 



2. If eq. (10251) is fulfilled, U is substituted by 



U'W = , ^ U'W. (C.29) 

VdetW 

Simulating the SU(2) Higgs model on a 12^ lattice with parameters /? = 2.0, 
K, = 0.25 and A = 0.5 (this point is in the confinement "phase") the change 



eq. (|C.29|) is accepted in 95% of the cases. 



C.2.2 Overrelaxation 

The new link U proposed is 

U' = WU^W = —^WU^W (C.30) 

This change is microcanonical, i.e. S{U) = — p/2tr (f/H^"!") = S{U') because of 
the reahty of the trace of SU(2) elements. Since the inverse relation of eq. ( |C.30| ) 
is U = WV^W , eq. ( |C.30| ) is a reflection of U and U' and therefore the factors pc 
cancel in the acceptance probability eq. ( p. 27 ). The change eq. ( C.30 ) is always 
accepted. 

C.3 Updating of the Higgs field 

For the updating of the Higgs field we use the algorithms proposed in |E3]. In 
particular the over-relaxation algorithm has been proved to be very efficient for 
the reduction of autocorrelations. 
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The part of the SU(2) Higgs model action eq. ( p.3| ) depending on the Higgs 
field reads 






n2 



(C.31^ 



Combining eq. ( |A.14| ) and eq. (|A.15|) the 2x2 matrix ip can be written in terms 
of the four real components 0i(x) {i = 1, 2, 3, 4) as 



(p{x) 



i)3{x) - i(j)i{x) 01 (x) + i(t)2[x) 
-01 (x) + 202(a;) 03(x) + i04(x) 



(C.32) 



If we do the same for the matrix 

B{x) = k'Y\U{x, fi)(p{x + jl) + U^{x — jl, fi)(p{x — p,) 



63 (x) — ib4^{x) bi{x) + ib2{x) 
-bi{x) + ib2{x) b3{x) + ib4{x) 



hen, (C.33) 

we can write the action eq. ( |C.31D in the four component notation of reference 



^l 



S^ = j:\i<l>i^)-bix)? + X 



(x)-l 



b\x) 



(C.34) 



where the squares mean the scalar product in R^. 

C.3.1 Heatbath 

We drop the index i and the argument x of the Higgs field variable 0i(x) {i = 
1,2,3,4) to be updated. A new Higgs variable 0' has to be generated according 
to the distribution 



dP(0') ~ e-^(<^') dV , ^(0) = (0 - bf + A(02 - 1)2 . 



(C.35) 



This is best achieved by introducing a free parameter a which parametrises the 
splitting of the potential V^(0) in a quadratic and in a quartic part 



\/(0) 



a(0-a-i6)^ + A(02-i;^)^-c„, 



where 



Cn 



1 + 



a — 1 



Xivt-l) + ia~'-l)b\ 



(C.36) 

(C.37) 
(C.38) 
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A trial Higgs variable is generated according to the Gaussian density (see eq. ( |C.6| )) 

2 



Ptriall0 j 

and is accepted with probability 



a 
- I e 

TT, 



■a(0'-a-lfe)2 



(C.39) 



-A(<^'2-t,2)2 



Wacptl</>J = e --- ^-' . (C.40) 

The parameter a is chosen so that it maximises the acceptance rate defined as 

.4(a) = y'dVptrial(0)^acpt(0), (C.41) 

which can be rewritten using eq. ( |C.36|) as 

A{a) = aV"M(l). 
Difi'erentiating with respect to a yields the cubic equation 



(C.42) 



/(a) = a^ ^ {1 - 2\)a^ - AXa 



2Xb' 



(C.43) 



The exact solution of eq. ( |C43| ) is inconvenient because the dependence on fe^ 
forces to solve it separately for each update. We use the approximate solution 
proposed in |^. The function /(a) in eq. ( |C.43| ) vanishes at 



tto 



2-^ + 



:i-A)^ + 4A 



1/2 



(C.44) 



Expanding /(a) in powers of {a — oq) and omitting the cubic term, eq. ( |C.43| ) 
becomes 



(ag + 4A)(a-ao) + (6ao + 4A-2)-(a-ao)^ = ^Xb"^ 



The positive solution of eq. ( 


C.45|) is 


a 


= /lo + [/il + /i2 &T^' 


with 




ho 


al + 4X 
"° 6«o + 4A - 2 


h 


/ al + 4X y 


h2 


4A 


6ao + 4A - 2 ■ 



(C.45) 

(C.46) 

(C.47) 
(C.48) 
(C.49) 
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The constants ao and hi (i = 0, 1, 2) can be computed in advance because they do 
not depend on h. For the update of it is only necessary to determine 6, find a 
according to eq. ( |C.46|) and the corresponding t;^ from eq. ( |C.37|) . Simulating the 
SU(2) Higgs model for relatively large value of A, the choice a = 1, corresponding 
to a naive splitting of the potential V{(j)) in eq. ( |C.36| ), would lead to a low 
acceptance A{1) in eq. ( |C.42| ). Using eq. ( |C.46| ) we get values A = 70% in the 
confinement "phase" (/5 = 2.0, k = 0.25, A = 0.5) and A = 59% in the Higgs 
"phase" (/9 = 2.3, n = 0.32, A = 1.0), both the simulations performed on a 12^ 
lattice. 



C.3.2 Overrelaxation 

The following change of the Higgs variable is proposed |^ 

— > (p' = 2a'^b - (p . 



(C.50) 



This change is not microcanonical, but it is a refiection of and 0', since the 
inverse relation of eq. ( |C.50|) is = 2a~^b — cj)'. According to eq. ( 2.27 ) the new 
Higgs variable is accepted with a probability 



PA 



min{l,exp [1^(0) -1/(0')]} 

min {l, exp [A(02 - 0'^) (0^ + 0'2 _ 2vl)] } 



(C.51) 



For the acceptance A of the change eq. (|C.50|) we get values A = 79% in the 
confinement "phase" (/? = 2.0, k = 0.25, A = 0.5) and A = 68% in the Higgs 
"phase" (/3 = 2.3, k = 0.32, A = 1.0), both the simulations performed on a 12"^ 
lattice. 
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Appendix D 

One-link integral and 
approximation of the modified 
Bessel functions 



In Sect. p.2| , we presented a method for reducing the statistical variance of correla- 
tion functions involving static charges which are represented in the path integral 
formalism by straight time-like Wilson lines. The time-like links can be replaced 
(with some restriction) by their expectation value (called one-link integral) in the 
fixed configuration of the other field variables. 



First of all, we present the derivation of the result eq. ( p.l9| ). We omit writing 



the (x, 0) dependence of U and W, for the rest the notation is the same as in 



Sect. 53 We denote by W the projection W/^Jdet{W) of W into SU(2). The 
one-link integral can be written after a change of integration variable as 

^ JdUUexp{{p/2)tT{U)} ^ 

/df/exp{(p/2)tr(f/)} ' ^ ■ ' 



where p = (3Jdet{W). To solve analytically these integrals, we use the following 
character expansion for SU(2) |^ 

exp{(p/2)tr(f/)} = -J2{2u+l)h.+,ip)x^^\U), (D.2) 

Pu 

where z/ = 0,1/2,1,3/2,... labels the irreducible representations of SU(2) with 
characters x^'^^ ^ind /„ is the modified Bessel function of integer order n. With 
help of eq. ( p. 2] ) and the orthogonality relations for the characters, it is easy to 
show that / dU exp{(p/2) tr (U)} = (2/p)/i(p). The integral in the numerator of 
eq. dDTD 

F = fdUUexp{{p/2)tT{U)} (D.3) 
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Figure D.l: Here, we show a histogram plot for the distrioution of p = 



pJdet(W{x, 0))) measured on each point x of a 6"^ lattice for 2000 field config- 
urations. The simulations are performed for two representative parameter sets. 



defines a linear transformation in C^ with the property that 

AF = FA, AgSU(2). 



(D.4) 



From the Schur's Lemma of representation theory, it follows that F is equal to 
a complex constant times the identity matrix 1. Writing the SU(2) link U in 
the quaternionic representation U = l/2tr(?7)l -\- iu ■ t, where Tj [i = 1,2,3) 
are the Pauli matrices and Ui G R, we see that only the component l/2tr(f/)l 



gives a contribution to F. By noting that tr (U) 



X 



(1/2) 



(f/) and using again 



eq. (p.2|) together with the orthogonahty property of the characters, we conclude 
that F = (2/p)/2(p)l, giving the desired result 



U 



Up) 
hip) 



w 



(D.5) 



The modified Bessel functions of integer order In{x) have the following asymp- 
totic behavior for large x 



iJx) 



V2 



-KX 



1+0 



X 



(D.6) 



For small x we get from the Taylor expansion of /„(x) 



In(x) 



1 /X 
d V2 



1 + {x') 



(D.7) 



In order to choose a good numerical approximation for the ratio R{p) = l2{p) / Ii{p) 
we measured the distribution of the argument p in Monte Carlo simulations on 



a 6 lattice for two values /3 = 2.0, 2.4. The results are shown in Fig. D.l. For 



p lying in the interval [0,13], we use two different Chebyshev approximations 
16^ , for the function R{p)/p (this choice is motivated by eq. (p.7|) ) in the range 



p G [0, 5] and for R{p) in the range p G [5, 13]. The Chebyshev approximation of 
a function f{x) in the interval x G [a, b] is defined as 



fix) - tc,Uy), y^^^i^^E [-1,1], (D.8) 

fc=0 2V" "') 



where Tk{y) = cos{k arccos(|/)) is the Chebyshev polynomial of degree k. The 
Chebyshev coefficients Ck, k = l,2,...,m are rapidly decreasing with k so that 
the error of the approximation is dominated by Cm+iTm+i{y) and is bounded by 
|cm+i|- We computed the coefficients Ck using the program MAPLE: the value 
of m is chosen such that we reach a precision of 10~^ in the considered ranges. 
We need m = 10 for the approximation of R{p)/p in the range p G [0,5] and 
m = 8 for the approximation of R{p) in the range p G [5, 13]. The Chebyshev 
polynomials satisfy the recurrence relation Tk+i{y) = 2yTk{y) — T^-iiy), k >1. 
An elegant and efficient way to evaluate the sum in eq. (p.8|) is the Clenshaw's 
recurrence formula described in reference [^ . 



When p > 13, which is rarely the case as can be seen from Fig. DJ^, we write 



Rip) ^ VPe-^W-2v/pe-/.(p)/P p_g) 

^e-^Hp) 

and use the polynomial approximations for ^/pe'^ In{p) {n = 0,1) in negative 
powers of t = p/3.75 given in reference ||100|| for the range p >= 3.75 (these 



approximations are motivated by eq. (p.6 )). For ^/pe '' Ii{p) we take the poly- 



nomial approximation up to the power t , which guarantees a precision of 10" 
for p > 13. Because of the extra power p~^ in the numerator of eq. ( p.9|) , we 
take the polynomial approximation for y/pe~^ Io{p) up to the power t~®, which 
also guarantees a precision of 10~®. 
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Appendix E 
Parallelisation 



In this appendix, we describe the parallehsation of the program for the simulation 
of the SU(2) Higgs model. Geometrically, the parallelisation consists in the parti- 
tion of the lattice on different processor elements (PEs). This is necessary for two 
reasons: the simulation is faster and the memory space needed on a 32*^ lattice 
to store the field variables does not fit into the 128 MB memory of one PE. We 
adopt a two-dimensional partitioning of the lattice. The xy-plane is distributed 
among the PEs and the t and z directions, for fixed x and y coordinates, are 
entirely contained on each PE. This choice was motivated by the measurement 
of the correlation function for the static potential: apart from smearing, it can 
be performed in the tz-plane on each PE independently. 

The programming aspects of the parallelisation are related to the communi- 
cation between PEs. The program is written in the language FORTRAN 90 for 
a CRAY T3E machine. For the communication between PEs we used the CRAY 
SHMEM (logically shared, distributed memory access) routines. In order that a 
local PE can pass (receive) a data object to (from) a remote PE there must be a 
known relationship between the local and remote address of the data object. This 
is realised by declaring the data object with the save attribute or in a common 
block. On CRAY systems, this means that the address of the data object in the 
memory of the local and remote PE is the same. These data objects are called 
symmetric. 

In this appendix, we set for simplicity a = 1. We emphasise all the names of 
variables, FORTRAN 90 statements and routines that we used in the program. 

E.l Geometry of the two-dimensional partition 



The lattice points in the xy-plane are partitioned among A^^: x Ny PEs. For fixed 
X and y coordinates, the points in t and z directions are entirely contained on 
the PE. Therefore, each PE works on a sub-lattice with L^ x Ly x L x T points, 

90 
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Figure E.l: Here, we show the xy-partition on each PE. The filled circles are the 
local points, i.e. the points on which the updating is performed. The empty circles 
are the remote points, i.e. the boundary points copied from the neighboring PEs. 
The dashed lines enclose points belonging to the same neighboring PE. The labels 
under the points are the values of the lexicographic index xs. 



where L^ = L/N^ and Ly = L/Ny. During the updating of the fields on a sub- 
lattice, we need field variables from sub-lattices having a common boundary in 
the xy-p\ane. Therefore, we store a copy of these boundary field variables. The 
number of points on each PE is then (L^ -|- 2) x [Ly + 2) x L x T. 

The PEs are labelled with a number npe that can be queried with the SHMEM 
routine my_pe{) and goes from to N^Ny — 1. To give an "identity" to the PEs 
(a relation with the geometrical sub-lattice they work on), we map the number 
npe into a "processor lattice" in the xy-plane by introducing the arrays 



peup{npe, i) and pedn{npe, 



(0 < npe < N^Ny - 1; i = 1, 2) , (E.l) 



where i = 1(2) labels the x{y) direction. The value of the array element peup{npe, i) 
{pedn{npe,i)) is the number of the nearest-neighbor PE in positive (negative) i 
direction of the PE labelled by npe. This is also the way periodic boundary 
conditions are implemented in the x and y directions. 



The xy-partition on a PE is represented in Fig. |E.1| . We decided to label the 
points on a sub-lattice with three coordinates 



{z, t, xs) , I < z < L , I <t <T , 1 < xs < volxy . 



(E.2) 



The index xs is the lexicographic index in the xy-plane and is constructed la- 
belling the local points (filled circles in Fig. [E.l| ), i.e. the points on which the 
updating is performed, from 1 to vol2 = L^Ly. The remote points (empty cir- 
cles), i.e. the boundary points copied from neighboring PEs, are labelled from 
V0I2 -|- 1 to volxy = {Lx + 2){Ly + 2). The order of the labelling is shown in 



Fig. |E.1|: the x coordinate runs fastest. 
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The space-time moveinents to the nearest neighbor points on the sub-lattice 
( "hopping" ) are given by the arrays 

zup{z) and zdn{z) , {I < z < L) , (E.3) 

tup{t) and tdn{t) , (1 < t < T) , (E.4) 

iup{xs,i) and idn{xs,i) , {1 < xs < volxy;i = 1,2) , (E.5) 

where i = 1(2) again labels the x{y) direction. The ending "up" {^^dn") refers to 
the positive (negative) direction. 

For each local point with index 1 < xs < vol2, we define an array boundary {xs) 
whose components are of the derived data type 

type, public :: boundaryjype 

logical :: tf 

integer :: ndir 

integer, dimension (3) :: xtg 

integer, dimension(3) :: pe 
end type boundary_type 

The component tf is .true, if the point xs is on the boundary of the x?/-partition 
(default value is .false.). The component ndir is the number of neighboring PEs 
that have a copy of the point xs: we have ndir = 3 for the points in the corners 
of the xy-partition and ndir = 1 for the other boundary points (default value is 
0). The component xtg{i) {i = 1, 2, 3) (default value 0) is the lexicographic index 
of the copy of the point xs on the neighboring (target) PE with number pe{i) 
(default value -1). 

The arrays (|E.1| ), ( [E.3|) , (|E.4| ) and ( |E.5|) together with boundary (xs) are de- 
fined in the module geometry. Once the geometry is initialised, it will be used 
throughout the program without modification. 

E.2 Communication between PEs 

The updating of the gauge and Higgs field requires the exchange of field variables 
between the PEs. This is the reason why we have a copy of the field variables 
at the boundary of the a;|/-partition belonging to neighboring PEs. It is essential 
that these boundary values are "up to date", i.e. they don't get modified simulta- 
neously on their local PEs. In order to avoid the use of field variables that do not 
have the current value, we associate to each PE a parity par which can assume 
the values or 1. According to the value of par, the updating is processed on 
points belonging to one of the two time sections 

T T 

l<t<- {par = 0) or - + l<t<T {par = 1) . (E.6) 

Before the updating starts, the value of par is initialised on the PEs. The PE 
with index has par = and the parity of two nearest-neighbor PEs is different. 
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Figure E.2: Here, we illustrate the updating of the gauge links. The filled circles 
represent points on which the links (full lines) are going to be updated. The empty 
circles represent points where the links (dashed lines) do not change because of 
the different parity of their PEs. The dashed boxes delimit the local points on the 
PEs. When PEO updates the links U{x, fi) with yU = 1 (/i = 2) it needs links with 
/i = 2 (yU = 1) from PE2 (PEl). If all PEs update links with the same direction 
fi and pass the links at the boundary of the a;|/-partition before changing the 
direction, there is no danger of using links from other PEs which are not up to 
date. 

Before changing the value of par, the PEs pass the field variables at the boundary 
of the xy-partition (in the time section which has been updated). This gives the 
correct algorithm for the Higgs field, which has only nearest-neighbor interac- 
tions. The updating of a gauge link on a point is more complicated since the 
plaquettes containing the link extend to next-to-nearest-neighbor points. When 
updating the gauge field in a given time section, the PEs must all work on the links 
U{x, n) with the direction n and they pass the boundary links before changing 
the direction. This procedure is necessary for two reasons. Firstly, when up- 
dating time-like links, gauge links from nearest-neighbor PEs at a different time 
coordinate are needed, but they are space-like. Secondly, the updating of hnks 
in the x{y) direction requires the use of gauge links at the same time coordinate 
on next-to- nearest neighbor PEs (which have also the same parity), but they are 
directed in the y{x) direction. This second situation is illustrated in Fig. [E.2| . 

The gauge field is represented by SU(2) matrices and the Higgs field can be 
represented according to eq. ( |A.15|) as a real constant times a SU(2) matrix. We 
can parametrise both of these types of 2 x 2 matrices as 




a = {a_re) + i{aJm), f3 = (bjre) + i{bJm) G 



(E.7) 



and we define for them the derived data type 
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type, public :: su2 

sequence 

real (kind= double) :: a-re,a_im,b_re,bJm 
end type su2 

The attribute sequence means that the components of the data object are stored 
in the sequence that they are declared. The name of the data type is su2 but it 
can be used for a larger class of matrices than only SU(2), as discussed above. 
The kind = double of the components means that they are double precision 
real numbers with a storage size of 8B. In the module su2_type we extend the 
intrinsic operators *,+,- to act between data objects of the type su2 with the 
usual meaning of matrix operations. We also extend the assignment = such that 
A = c sets the components of a su2 data object A corresponding to the identity 
matrix times the real number c. We define new operators between su2 data 
objects .doth, and .hdot., where A.doth.B and A.hdot.B correspond to the matrix 
multiplications AB"^ and A^B respectively. Moreover, we define the function det, 
which has a su2 data object as argument and returns the value of the determinant 
of the corresponding matrix. There is one point concerning the precedence of 
the operators which is worthwhile to mention. An extended intrinsic operator 
maintains its precedence, a defined binary operator has the lowest precedence. 
For example, the expressions 

A.doth.B + C and {A.doth.B) + C (E.8) 

correspond to the matrix operations A{B + C)'^ and AB'^ + C respectively, which 
are clearly different. 

The gauge and Higgs field on the lattice are represented by arrays u and phi 
of data objects of the type su2 with indices 

u{z,t,xs,niu) and phi{z,t,xs) . (E.9) 

The ordering of the indices is chosen carefully. The language FORTRAN stores 
the array A{ii,i2, ...,in) in the computer memory by columns. This means that 
the index ii is the fastest and the index in the slowest. When the array elements 
are used in the program they are copied into the memory cash of the machine. 
A better performance of the code is reached when the sequence of the indices 
ii, ...,in in the definition of the array A corresponds to the inverse sequence of do 
loops over the indices in the program (the loop over in is the outermost and the 
loop over ii the innermost). Why the sequence of indices in ( |E.9| ) is the best for 
our purposes will become clear soon. 

The communication between PEs on a CRAY T3E machine is implemented 
with the SHMEM routines shmemjput and shmem_get. The routine shmemjput 
has the following arguments: 

shmem_put (target, source, size, target PE j . 
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The first argument target is the target array on a remote PE to which we put 
the data. Actually, what is passed in FORTRAN is the address in memory of the 
first element of the array target. The array target must be a symmetric array, 
i.e. its components have the same address in memory on each PE. The second 
argument source is the array on the local PE from which the data are transferred. 
The third argument size is the number of words w (1 w = 8B = 64 bits) that 
are transferred. The last argument is the number of the remote PE to which the 
data are passed. The arguments of the routine shmem_get are 

shmem_get (target, source, size, source VE) . 

The interpretation is similar but the target array, to which the data are trans- 
ferred, resides on the local PE and the source array, from which we get the data, 
is on a remote PE, whose number is given by the last argument. The array source 
must be symmetric. 

As concerns the time duration of the communication between PEs, a useful 
quantity to know is the time interval between the start of the data transfer and 
the availability of the data on the target PE, called data latency tl. For the 
CRAY T3E we have 

Tl = O.SxlO-^'s. (E.IO) 

In a simplified but rather realistic model, we can parametrise the total time per 
byte needed for transferring a block of data of size P (in bytes B) as 

tB = a + -, (E.ll) 

where a is the asymptotic (P -^ cxd) transfer rate (bandwidth) and has the value 
1 s/300 MB. Identifying the data latency as the time for transferring 1 B, we can 
set the parameter h k, n^. If we define the critical size P^ for the block of data to 
be transferred from the condition 

ha , , 

— = — , E.12 

Pc 10 ^ ^ 

which means a 10% deviation from the asymptotic rate of transfer, we get 
Pc = 1000 B. For this size of the data block transferred, the time losses due 
to communication are only 10%. 

The fundamental operation of synchronisation of the PEs is achieved by calling 
the SHMEM routine 

barrier (). 

The call of the barrier {) routine "announces" the arrival of a PE at that line 
of the code. The PE suspends the execution of commands until all PEs have 
called barrier {). A barrier ensures that a PE, prior to synchronising with other 
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PEs, has completed all previously issued local memory stores and remote memory 
updates issued via SHMEM routine calls such as shmenijput. 

In the following example of parallelised code, we give the structure of the 
updating of the gauge field contained in the module update^gauge field. The 
outermost loop changes the parity par of the PE which in turn determines 
through ( |E.6| ) the time section tmin{par) < t < tmax{par) in which the up- 
date is performed. The second loop is over the directions mu of the links, for 
the reasons we explained at the beginning of this section and also because of 



the integrated autocorrelation times as discussed in Sect. |2]J. The third loop 
runs over the lexicographic index xs in the local x^-plane. For fixed par, mu 
and xs, the links u{z,t,xs,mu) in the specified range of t and for all z are up- 
dated: if xs is at the boundary of the a;?/-partition, the package of field variables 
■^(1 < z < L, tmin{par) < t < tmax{par),xs, mu) is sent to the neighboring PEs 
that need a copy of it. The information on which PE needs a copy and where 
the copy has to be stored is contained in the components of the array element 
boundary{xs). The parameter package is the size in words of the data block to 
be sent and is equal to 4 x L x T/2. For L = 32 = T this corresponds to 16 KB, 
which is well above the critical size Pc = 1000 B in eq. (Pi).12| ). Now, it becomes 
clear why we choose the order given in ( |E.9| ) for the arguments of the fields: the 
z coordinate is always sent entirely in the package, the t coordinate between a 
lower and an upper bound. This package must lie in one memory sequence, in 
fact what is passed to the target PE is the address of the first component of the 
sequence (as can be seen in the code example) and the length in words of the 
sequence. The argument m,u of the gauge field in (|E.9| ) is the last because of the 
structure of the updating. After having updated the links in a time section for 
fixed direction mu, there is a barrier. The copies of the boundary links must be 
up to date before processing the next direction. Once the updating in the time 
section is terminated, the parity of the PE is changed. The updating of the Higgs 
field in the module updateJiiggs field has the same structure without loop over 
m,u. 

Actually, not all the boundary field variables that are copied are needed during 
the updating. The simplified structure and the relatively large size of data blocks 
transferred justify these small losses in performance. The performance of the 
whole program has been established with the tool APPRENTICE to be 

93.4Mfiop/s/PE. (E.13) 

For comparison, the theoretical performance of the CRAY T3E machine is 900 Mfiop/s/PE 
but ( PH.ISD is a respectable number for a program written in a high level language. 
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Example of parallelised code: updating of the gauge field. 

do pp=l,2 

do mu=0,3 

do xs=l,vol2 

do t=tmin(par) ,tmax(par) 
do z=l,L 

! updating of u(z,t,xs,mu) 
end do 
end do 

! transfer of fields on boundary points 
if (boundary (xs) %tf. eqv. . true.) then 
do i=l, boundary (xs)%ndir 

call shmem_put(u(l,tmin(par), boundary (xs)%xtg(i),mu), 



& 



end do 
end if 

end do 

! synchronisation 

call barrier() 

end do 

! parity switch 

par=modulo (par+ 1 , 2) 

end do 



u(l,tmin(par),xs,mu), & 
package, & 
boundary (xs) %pe (i)) 



E.3 Program overview 

Our program for the simulation of the SU(2) Higgs model on the lattice is com- 
posed of the following modules: 
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module su2_type 

module global 

m,odule geom,etry 

module ranfloat 

module initiaLconfig 

m,odule run_updating 

module update^gaugefield 

m,odule update Jiiggsfield 

m,odule sm,earing 

module observables 

m,odule kinds 

module timing 

module utilities 

module gauge-transformation 

program main 

Some of the modules have been already described. Here, we give a short descrip- 
tion of the other modules in which important definitions or subroutines of the 
program can be found. We refer to the comments in the code for a more detailed 
information. 

The module global contains the global variables and parameters. The global 
parameters for which the user can choose a value (with some restrictions) are: 
the numbers of PEs in x and y directions (A^^ = npe_x and A^^ = npejy: must 
be either even numbers or 1), the lattice sizes L = spacesize (must be divisible 
by Nrc and Ny) and T = time_size (must be even), the numbers of smeared 
fields {smleveljphi and smlevelJink for the Higgs and gauge field respectively) 
and corresponding smearing levels (arrays philevels and ulevels), the smearing 
strength e = omegaJink for the gauge field (see eq. (|4.4| ), the smearing procedure 
for the Higgs field in eq. ( p.23| ) has no free parameters) and the maximal space and 
time extensions of the correlations for the measurement of the static potentials 
and the static-light meson spectrum {rmax and tmax). 

The module ranfloat contains the random number generator as proposed by 
M. Liischer PB| , ^^ and subroutines for the generation of random numbers with 



distributions needed in the updating algorithms, see Sect. |C.l 



The module run_updating organises the updating of the fields, the measure- 
ments of the observables and writes the measurements into files. 

The module smearing contains the computation of the one-link integrals (see 
Appendix 0) and the construction of the smeared fields. The subroutine ini-eigphi{) 
contains the coefficients describing the approximate wave functions for the static- 
light meson ground, first and second excited states. 

The module observables contains the subroutines for the measurements of the 
plaquette, the Higgs length squared and its square, the gauge invariant links, the 
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matrix correlation functions for the static potentials (in the tz-plane) and for the 
static-light meson spectrum. The translation invariance and the isotropy of the 
space are used to average the measured quantities over the lattice points and the 
directions. 

In order that the program main is executed, several run parameters are read 
during execution from a file parameter, d. The run parameters needed are listed 
at the end of this section. For the continuation of an old run {start = 1) the 
parameters from 4. to 9. must be omitted, for the start from a given configuration 
{start = 2) they must be substituted by 4. confpath (name of the path where 
the start configuration is stored) and 5. confname (basic name of the start 
configuration files). The parameter bitseed is the initialisation of a random bit 
generator |^ which produces for each PE different initialisation seeds of the 
Marsaglia-Zaman random number generator. The random number sequences on 
the PEs are different for different seeds IPBI. We recall that one iteration is 



composed of one heatbath sweep for both gauge and Higgs field and Nqr times 
the combination of one over-relaxation sweep for the gauge field plus three over- 
relaxation sweeps for the Higgs field. The measurements of the matrix correlation 
functions for the static potentials and for the static-light meson spectrum are 
performed only each iterjtojneas iterations. 
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Parameters to be written in the file parameter. d. 

1. path 

full path where to write the output files 

2. name 

basic name for the output files 

3. start=0,l,2 

0: new run, 1: continuation of old run, 2: start from a thermalised 

configuration 

Ji^. P (or confpath) 

5. K, for confname) 

6. A 

7. 1 < hitseed < 2^3 - 1 

seed for the initialisation of the random number generators 

8. order 

0: hot start, ^ cold start 

9. therm 

number of thermalisation iterations 

10. NJter 

number of updating iterations 

11. cpuMme-max 
maximal CPU time in seconds 

12. N_OR 

ratio (over-relaxation sweeps) /(heatbath sweeps) for the gauge field 

13. iter^to-meas 

ratio (iterations up dating) /measurements 

14- N_save 

number of iterations after which the data are written into files 
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